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Abstract 


An  analytical  study  of  a  technique  for  achieving  high 
Reynolds  number  similarity  in  a  wind  tunnel  was  conducted. 
The  technique  under  consideration  required  that  the  momen¬ 
tum  boundary  layer  be  immersed  in  a  layer  of  eryogenically 
cooled  gas  so  that  boundary  layer  development  would  be 
determined  by  the  low  temperature  fluid  properties.  The 
experimental  technique  was  shown  to  be  theoretically  feas¬ 
ible,  with  effective  Reynolds  number  increased  by  as  much 
as  a  factor  of  seven  or  eight,  provided  that  sufficiently 
high  mass  rates  of  flow  of  cool  gas  could  be  provided. 


AN  ANALYTICAL  STUDY  OF  A 


LOCALLY  COOLED  HIGH  REYNOLDS  NUMBER 
SIMULATION  TECHNIQUE 


I.  Introduction 


Aerodynamic  structures  are  often  tested  in  wind  tunnels 
to  predict  their  free  flight  performance.  The  aerodynamic 
forces  on  a  structure  strongly  depend  on  the  Reynolds  number, 
so  it  is  important  that  the  free  flight  value  of  this  para¬ 
meter  be  matched  as  closely  as  possible  by  tunnel  test  con¬ 
ditions.  Due  to  the  limited  scale  size  of  models  that  are 
used  in  wind  tunnel  testing,  however,  this  requirement  can¬ 
not  always  be  properly  met.  Many  free  flight  conditions  of 
interest  have  not  been  adequately  simulated  using  current 
techniques . 

There  is  a  desire  to  develop  new  procedures  which 
might  remove  some  of  the  existing  constraints  on  wind  tun¬ 
nel  Reynolds  numbers.  The  study  of  this  paper  is  concerned 
specifically  with  a  proposal  to  simulate  high  Reynolds 
number  flight  conditions  by  locally  immersing  the  model 
being  tested  in  an  envelope  of  cold  fluid.  Because  the 
Reynolds  number  increases  significantly  as  temperature 
decreases,  it  is  hoped  that  such  a  technique  could  increase 
the  effective  wind  tunnel  Reynolds  number  by  as  much  as 
an  order  of  magnitude. 


Background 


Wind  tunnel  testing  is  basea  on  the  premise  that 
the  fluid  flow  around  a  scale  model  of  a  structure  will  be 
similar  in  some  sense  to  a  fluid  flow  around  the  full  scale 
structure  itself.  In  the  case  where  the  fluid  flow  is 
inviscid  and  incompressible,  this  premise  would  allow  the 
results  of  a  tunnel  test  to  be  extended  to  an  arbitrary 
situation  as  a  function  of  the  scale  of  the  model  and  of 
the  free  stream  velocity  of  the  fluid. 

Unfortunately,  real  fluids  are  neither  inviscid  nor 
incompressible.  This  is  a  significant  comnlicat ion ,  since 
compressibility  and  viscosity  are  dependent  on  the  partic¬ 
ular  conditions  for  which  they  are  being  measured.  Con¬ 
sequently,  it  is  not  immediately  apparent  how  results  based 
on  the  effects  of  viscosity  and  comppressibility  in  a  wind 
tunnel  can  be  used  to  make  predictions  about  the  behavior 
of  a  structure  in  the  general  case. 

In  order  to  gain  insight  into  the  question  of  how 
wind  tunnel  test  results  can  be  extended  to  other  situa¬ 
tions,  it  is  convenient  to  non-dimensionalize  the  equations 
which  describe  the  behavior  of  real  fluids.  When  this  is 
done,  two  parameters  appear  to  be  of  major  significance 
when  the  fluid  flows  around  two  geometrically  similar 
structures  are  to  be  compared. 

The  first  of  these  parameters  is  the  Reynolds  number. 
Denoted  by  Re,  and  determined  from  the  relationship  Re  = 


r.  .’.v  v-  .v  .  .  -v  .r 


W^i/1  Jr| 


^  the  Reynolds  number  is  a  measure  of  the  ratio  between 
the  inertia  force  of  the  flow  and  the  force  which  results 


from  viscous  effects. 

The  second  parameter  is  the  Mach  number,  which  is 
determined  from  the  relationship  M  =  ^.  Commonly  known 
as  the  ratio  of  local  velocity  to  local  speed  of  sound,  the 
Mach  number  is  a  measure  of  the  way  in  which  a  disturbance 
in  the  flow  field  (an  airfoil,  for  instance)  may  effect  the 
fluid  flow  upstream  from  the  disturbance. 

When  the  aerodynamic  forces  on  a  structure  are  non- 
dimensionalized  with  respect  to  dynamic  pressure  (q  =  <spV2) 
and  to  some  characteristic  area,  theoretical  and  empirical 
results  suggest  that  these  non-dimensionalized  forces  are 
functions  primarily  of  the  Reynolds  and  Mach  numbers.  There 
fore,  the  aerodynamic  forces  on  a  structure  can  be  used  to 
predict  the  aerodynamic  forces  on  a  geometrically  similar 
structure  if  the  Mach  and  Reynolds  numbers  for  the  two  flow 
situations  are  the  same. 

In  principle,  the  behavior  of  an  aerodynamic  struc¬ 
ture  in  free  flight  can  be  predicted  from  wind  tunnel  tests 
on  a  scale  model  if  the  tunnel  Mach  number  and  Reynolds 
number  can  be  set  equal  to  the  free  stream  values  of  those 
parameters . 

This  principle  cannot  often  be  applied  in  practice. 

Wind  tunnels  are  limited  in  the  range  of  Mach  and  Reynolds 
numbers  over  which  tests  can  be  conducted,  and  it  is  rarely 
possible  to  simulate  the  free  stream  values  of  these  two 


parameters  simultaneously. 


There  are  some  flow  situations  of  practical  interest 
where  the  dependence  of  aerodynamic  forces  on  the  Mach 
number  is  negligible  (when  M  =  0.3-0. 5  or  less,  see  ref  1). 
In  these  cases,  wind  tunnel  test  results  can  be  used  to  give 
good  predictions  of  free  flight  behavior  if  Reynolds  number 
similarity  can  be  obtained  in  the  wind  tunnel. 

Even  in  this  case  where  Mach  similarity  is  insignif¬ 
icant  there  are  substantial  constraints  on  wind  tunnel 
testing . 

In  an  attempt  to  simulate  high  Reynolds  number  con¬ 
ditions,  the  wind  tunnel  experimentalist  may  try  to  adjust 
the  free  stream  tunnel  velocity  (V),  the  characteristic 
length  of  the  model  (l) ,  or  the  kinematic  viscosity  of 
the  fluid  ( v  =  u/p ) . 

The  tunnel  velocity  is  limited  by  the  ability  of  the 
tunnel  to  accelerate  the  flow.  There  is  an  upper  limit 
to  the  velocity  that  can  be  achieved  in  the  test  section  of 
a  wind  tunnel.  Furthermore,  high  tunnel  velocities  result 
in  high  test  Mach  numbers,  so  that  compressibility  effects 
have  a  significant  impact  on  the  test  results  in  spite  of 
the  fact  that  compressibility  may  not  be  important  in  the 
free  fligh  situation  to  be  simulated. 

The  characteristic  length  of  the  model  is  obviously 
limited  by  the  size  of  the  tunnel  test  section.  Although 
large  tunnels  could  be  built,  tremendous  amounts  of  power 
would  be  required  to  drive  such  a  tunnel. 
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The  kinematic  viscosity  has  a  rather  strong  dependence 
on  pressure  and  temperature.  Kinematic  viscosity  decreases 
with  increased  pressure  or  decreased  temperature,  so  high 
pressures  or  low  temperatures  can  be  used  to  obtain  high 
Reynolds  number  simulation  in  a  wind  tunnel. 

The  National  Transonic  Facility  (NTF) ,  currently  being 
constructed  for  NASA  at  Langley,  Virginia,  will  combine 
increased  pressure  and  cryogenic  temperatures  in  order  to 
simulate  high  Reynolds  number  flow  (ref  2) . 

There  are  other  techniques  which  can  be  used  to  sim¬ 
ulate  high  Reynolds  number  flow.  Since  such  flows  are 
generally  characterized  by  turbulent  boundary  layers,  tur¬ 
bulence  can  be  induced  over  the  structure  to  be  tested 
(through  use  of  roughness  or  "trip  wires,"  for  example;  see 
ref  3) ,  thereby  simulating  high  Reynolds  number  flow  in  an 
otherwise  low  Reynolds  number  situation.  Another  tech¬ 
nique  involves  testing  several  different  size  models  with 
similar  geometries  and  attempting  to  extrapolate  the 
results  from  these  tests  to  appropriately  sized  structures. 

Each  of  the  techniques  described  above  has  drawbacks, 
typically  with  respect  to  added  expense  of  the  test  proced¬ 
ure,  added  complexity  of  the  scale  model  to  be  tested,  or 
theoretical  difficulty  in  extending  test  results  to  the 
free  flight  situation  of  interest.  Therefore,  a  relatively 
simple,  low  cost  means  of  providing  high  Reynolds  number 
simulation  is  clearly  desirable. 

On  possibility  for  such  a  technique  has  been  sug- 


gested  by  W.  Luchuk  of  the  Air  Force  Arnold  Engineering 
Development  Center  (ref  4) •  Mr.  Luchuk  has  proposed  that 
high  Reynolds  number  simulation  could  be  obtained  locally  by 
injecting  a  cryogenic  fluid  to  envelop  the  momentum  boundary 
layer . 


Purpose 

The  purpose  of  this  study  is  to  analyze  the  potential 
of  a  high  Reynolds  number  simulation  technique  which  accom¬ 
plishes  boundary  layer  cooling  through  the  introduction  of 
a  cool  gas  over  the  structure  to  be  tested. 

In  this  paper,  the  major  considerations  that  must  go 
into  the  development  of  a  locally  high  Reynolds  number 
simulation  technique  are  outlined  and  discussed.  These 
considerations  are  then  applied  to  establish  the  important 
physical  parameters  for  such  a  technique  and  to  make  con¬ 
clusions  regarding  the  feasibility  of  the  technique. 

The  flow  situation  of  interest  is  illustrated  in 
fig  1,  In  particular,  the  case  for  flow  over  a  flat 
plate  will  be  studied  in  detail. 

The  results  of  the  analysis  described  above  are 


discussed  at  the  end  of  this  paper,  and  recommendations 
for  additional  study  are  also  included  there. 


II .  Theor 


The  mathematical  model  which  is  used  to  describe  the 
behavior  of  a  fluid  is  typically  developed  from  an  Eulerian, 
or  "Control  Volume,"  point  of  view.  From  this  point  of  view, 
the  laws  of  physics  which  require  conservation  of  mass,  mornen 
turn,  and  energy  of  a  fixed  mass  system  are  adapted  to  predict 
behavior  in  and  around  a  fixed  region  in  space. 

The  Eulerian  analogs  to  the  fixed  mass  concepts  of 
conservation  of  mass,  Newton's  Second  Law,  and  conservation 
of  energy  are,  respectively,  the  Continuity  Equation,  the 
Momentum  Equation,  and  the  Energy  Equation. 

These  equations  are  presented  below  in  two  forms. 

In  the  first  case,  they  are  presented  in  tensorial  form 
so  that  they  may  be  applied  to  any  flow  situation  describ- 
able  in  terms  of  an  orthogonal  coordinate  system.  The 
principles  which  are  applied  to  the  case  of  a  flat  plate 
can  then  be  extended  to  apply  to  arbitrary  shapes  when 
an  appropriate  coordinate  system  has  been  defined. 

In  the  second  case,  the  equations  are  presented  in  the 
well  known  form  for  rectangular  (Cartesian)  coordinate 
systems . 

Appendix  A  to  this  report  gives  a  detailed  mathematical 
derivation  of  the  governing  equations  from  a  tensor  analysis 
point  of  view. 
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The  Continuity  Equation 


The  amount  of  mass  contained  within  a  control  volume 
may  change  with  time.  This  change  can  be  the  result  of  a 
change  in  the  density  of  the  mass  contained  within  the 
volume  due  to  an  influx  of  mass  to  the  control  volume.  This 
is  expressed  as 


+  p(V-v)  =  0 


In  terms  of  Cartesian  coordinates,  with  v  =  ui  +  vj  +  wK, 
this  is  written 


l£  +  ul£  +  vl£  +  wj£  +  D(lH  +  IX  +  lX)  =  o  (2) 

3t  +  u3x  V3y  W3  z  pk3x  3y  3z'  U 


In  the  case  where  density  can  be  assumed  to  be 
constant,  these  equations  reduce  to 


(V»v)  = 

3u  ,  3v  ,  3w 
3x  3y  3z 


The  Momentum  Equation 


The  momentum  within  a  control  volume  can  change  with 
time  due  to  changes  in  properties  in  the  control  volume 
caused  by  the  transfer  of  momentum  into  the  control  volume 


'  *  1  *  *j>  “ji  *  '  •  *  *  k'*  A  ,*•  ,* 


9 


or  due  to  the  action  of  forces  on  the  mass  inside  the  volume 


In  the  case  where  the  shear  stresses  acting  on  the 
boundary  of  the  control  volume  are  linearly  dependent  on  the 
rate  of  deformation  of  the  mass  contained  within  the  volume, 
the  momentum  equation  is  given  by 


p-^  =  -Vp  +  V(XV»v)  +  V»y(Vv  +  grad  v)  +  pf 


(5) 


Equation  (5)  is  really  a  system  of  three  scalar 
equations,  the  Navier-Stokes  equations,  which  are  written 
in  Cartesian  coordinates  as 


x-direction: 
3u 


D(-°“  +  u—  +  v—  +  w— )  =  X  -  & 
PV3t  u3x  V3y  3z;  A  3x 


3  r  ( o  3u 

"  3^2— 


3x 


%  (l¥  +  f  +  |H))} 


3x 


3u> 


y-direct ion : 


/ 3v  .  3v  ,  3v  ,  3v  \  „ 

p(n  +  +  ^  +  v~z]  =  1 


l£ 

9y 


,  d  f  ( o  dV  2/  (  i  ±  ^  V  i  dW\  n  I 

+  -  /a  (3^  +  37  +  3i}^ 


,  o  r  /  dV  ,  d q  ,  d  r  /  dV  ,  dU\i 

+  +  +  s^<u(^  +  w’1 


(6) 


z-direction : 


In  the  case  where  body  forces  can  be  neglected  and 
properties  can  be  assumed  constant,  these  equations  become 


=  -Vp  +  pV2v 


x-direction: 


/ 3u  .  3u  ,  3u  .  3u\ 

p(Tt  +  +  +  wsi>  = 


y-direct ion ! 


r  3v  ,  3v  ,  3v  ,  3v\ 

p  (tt  +  U-T—  +  V-r —  +  w-r— ) 
3t  3x  3y  3z 


z-direction: 


/  3w  ,  3w  ,  3w  ,  3w\ 

p(n  +  +  +  "sf  = 


3p  ,  1 32 u  ,  9 2 u  ,  32 
-at  +  y(3^  +  3^  +  3i 


3p  ,  / 32 v  ,  32 v  ,  32 

■ay  +  w(a^  +  3^  + 


3p  ,  /32w  ,  32w  ,  32 

■aS  +  w(T3F  +  3^  +  15 


The  Em 


Equation 


The  energy  contained  within  a  fixed  region  of  space 
may  change  with  time  due  to  the  addition  of  mass  or  heat 
to  the  control  volume  and  due  to  work  done  on  the  control 
volume.  When  this  statement  is  expressed  mathematically, 
the  resulting  equation  can  be  combined  with  the  momentum 
equation  to  yield  the  thermal  energy  equation 


Pfrr  =  -pV*v  -  V*q  +  $ 


For  a  perfect  gas,  and  when  viscous  dissipation  is 
negligible,  equation  (9)  reduces  to  the  form 


(10) 


DT  Dp  _  - 
pCpDt  =  Dt  '  7*q 


Application  of  Fourier’s  Law  and  the  assumption  of 
constant  properties  gives  the  familiar  form 


T5S  '  “v2T 


which  is  written  in  terms  of  Cartesian  coordinates  as 


3T  ,  3T  ,  3T  ,  3T  ,32T  ,  32T  ,  32T\  ,10s 

at  +  ual  +  v37  +  ”ai  =  “<3^  +  aF  (12) 
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Reynolds  Similarity  and  Wind  Tunnel  Testin 


From  the  presentation  above  it  is  clear  that  the 
momentum  equation  is  coupled  to  the  energy  equation  since 
the  fluid  properties  are  dependent  on  the  thermodynamic 
state  of  the  fluid.  In  terms  of  its  effect  on  wind  tunnel 
testing,  it  is  this  coupling  which  results  in  the  require¬ 
ment  that  tunnel  values  of  Re  and  M  should  be  simultaneously 
set  equal  to  their  free  flight  values. 

In  many  practical  applications  it  is  reasonable  to 
treat  the  fluid  properties  as  essentially  constant.  In 
this  case  the  momentum  equation  is  uncoupled  from  the 
energy  equation  and  the  velocity  profile  of  the  fluid  can 
be  obtained  from  the  momentum  equation  alone.  Since 
aerodynamic  forces  on  a  structure  are  dependent  on  the 
velocity  profile  of  the  fluid  surrounding  it,  the  aero- 


H  -rm  q  m  t  o  o  ut  1  1  "K  o  n  D-nonrlorrl'  -p  Vi  o  rvrnnDT+.'i  oc  P 


These  non-dimensional  terms  are  denoted  with  a  super¬ 
script  cross  (+)  and  are  defined  as 


t  =  t /U/MJ 

v*  =  v/V. 

P+  =  P/P„ 


p+  =  (p-p  )/p  v2 

1  ~  rro  ^  nn  on 

+ 


x  .  =  x  .  / 1 

co  oo  oo  j 


A'  -  a/ A 
'+  =  fV 


u  =  u/u  (13) 


Substituting  for  the  dimensional  quantities  in 
jquation  (5),  the  result  is 


— + 


p  V2  , Dv 

oo  oo  + 

l  Dt" 


p  V2  A  V  ...  ,  U  v 

00  oo  +  +  CO  oo  -f  .  - T  \  oo  oo  -f  .  -f  4-  —  4-  . 

- v  p  +  - V  ( A  V  •  v  )  +  - V  -(u  def  V  ) 


(14) 


where  def+v+  =  V  +  v  +  +  grad+v+. 

Combining  terms  and  simplifying  gives  the  result 


P  — -r-  =  -V  +  p+  +  (1/Re  )  V+(  A+V+*v+)  +  (1/Re  )V+*(p  +  def  +  v+) 

Dt  ,A 


V  p  l  V  p  l 

i  i-.  OO  OO  Oo'OT 

where  Re  ,  =  — ; -  and  Re  =  - 

co ,  A  A  °o  p 


(15) 


The  quantity  Re  is  the  Reynolds  number. 

In  the  case  of  constant  properties,  p+,  y+ ,  and  A+ 
are  all  equal  to  unity,  so  the  non-dimensional  momentum 
equation  takes  the  form 


Dt 


~~-r  =  ( 1  /Re  -,)V+(V+*v+)  +  ( 1  /Re  )  V+  •  (def  +  v+  )  -  V+t>  + 


"oo  ,  \  ‘ 


(16) 


If  two  flows  have  the  same  geometry  (so  that  the 
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pressure  gradient  and  boundary  conditions  are  the  same  for 


both)  and  if  the  assumption  of  constant  properties  applies 
equally  to  the  two,  they  will  have  the  same  solution  to 
equation  (16)  if  they  have  the  same  value  of  Reynolds  number 

When  a  fluid  in  motion  encounters  a  solid  surface,  the 
fluid  and  the  surface  will  exert  forces  on  each  other.  The 
magnitude  of  these  forces  is  related  to  the  development  of 
the  boundary  layer  (the  region  in  which  the  solid  object 
has  a  significant  effect  on  the  flow  field).  It  is  clear 
from  equation  (16)  that  the  development  of  the  boundary 
layer  is  dependent  on  the  Reynolds  number,  so  the  shear 
forces  a^e  dependent  on  Re  as  well.  In  particular,  when 
the  shear  stress,  x,  is  non-dimensionali zed  in  terms  of 
a  friction  coefficient,  c„  =  x/(gp  V2  )  ,  this  coefficient  is 
a  function  of  the  Reynolds  number  (ref  1). 

If  two  flows  are  geometrically  similar  and  if  they 
have  the  same  Reynolds  number,  they  will  have  the  same 
value  of  c^. 

Wind  tunnel  testing  is  based  on  this  concept  of  sim¬ 
ilarity.  A  scale  model  of  a  structure  can  be  placed  in 
a  wind  tunnel  and  values  of  c^  as  a  function  of  Re  can 
be  determined.  These  values  of  c^.  can  then  be  applied  to 
predict  the  behavior  of  the  full  scale  structure. 

Unfortunately,  the  concept  of  Reynolds  similarity 
breaks  down  for  high  values  of  Re  and  simulation  of  high 
Reynolds  numbers  typically  is  not  achieved.  Chapter  I 
of  this  report  has  detailed  some  of  the  reasons  for  the 


failure  of  this  similarity  concept  and  has  suggested  the 
desirability  of  designing  techniques  to  permit  high  Rey¬ 
nolds  number  simulation. 

To  obtain  high  Reynolds  numbers  on  sub-scale  models, 
attention  has  been  focused  on  the  possibility  of  varying 
the  quantity  v  (v=  y/p)  which  appears  in  the  definition 
of  Re.  Since  the  value  of  v  decreases  with  a  decrease  in 
temperature  for  a  gas,  cooling  techniques  can  be  used  to 
provide  high  Reynolds  number  simulation. 

The  remainder  of  this  paper  is  concerned  with  one 
possible  technique  for  the  use  of  cooling. 

If  a  cool  gas  can  be  introduced  into  the  tunnel  free 
stream  so  that  the  structure  being  tested  is  enveloped 
entirely  by  the  cool  fluid,  the  Reynolds  number  in  the 
vicinity  of  the  structure  will  be  significantly  greater 
than  the  Reynolds  number  based  on  properties  of  the  tunnel 
free  stream  fluid. 

Results  to  be  presented  elsewhere  in  this  paper 
will  suggest  that  the  use  of  sufficiently  cooled  fluids 
could  increase  the  effective  Reynolds  number  by  a  factor 
of  seven  or  eight. 

The  feasibility  of  using  this  technique  of  locally 
high  Reynolds  number  simulation  will  be  discussed  and 
analyzed  in  the  paragraphs  that  follow. 


The  high  Reynolds  number  simulation  technique  which  is 
the  subject  of  this  study  requires  the  immersion  of  the 
structure  being  tested  in  an  envelope  of  cool  fluid.  If 
the  momentum  boundary  layer  develops  entirely  within  the 
cool  fluid  layer,  the  Reynolds  number  which  determines  the 
non-dimensional  force  coefficient  will  depend  on  the  prop¬ 
erties  of  the  cool  fluid. 

In  particular,  the  ratio  of  Reynolds  number  based  on 

cool  fluid  properties  to  Reynolds  number  based  on  tunnel 

fluid  properties  can  be  estimated  by  Re  /Re  ~  v  /v  . 

c  s  sc 

For  cryogenic  temperatures,  the  ratio  v  /v  may  be  as 

s  c 

high  as  seven  or  eight  (see  Appendix  D) . 

In  order  to  take  advantage  of  the  potential  for  such 
an  increase  in  effective  Reynolds  number  it  is  necessary  to 
study  the  ways  in  which  the  tunnel  fluid  may  interact  with 
the  cool  fluid. 

Three  major  fluid  interactions  most  significantly 
impact  the  development  of  a  Reynolds  number  simulation  tech 
nique  involving  the  injection  of  a  cooled  gas.  These 
are  momentum  mixing,  thermal  mixing,  and  condensation/icing 
Momentum  mixing  will  occur  when  fluids  moving  paral¬ 
lel  to  each  other  have  different  velocities  (see  fig  2). 

The  shear  stress  created  at  the  velocity  discontinuity  will 
act  to  smooth  the  discontinuity  and  will  result  in  the 
development  of  a  momentum  mixing  region. 


v\ «. 
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Figure  3.  The  Thermal  Mixing  Region 


v 

V 


Thermal  mixing  occurs  when  a  fluid  in  motion  is 
separated  into  regions  held  at  different  temperatures  (see 


fig  3) .  The  temperature  discontinuity  will  create  a  thermal 
gradient  and  the  consequent  development  of  a  thermal  mixing 
region . 

In  the  case  where  a  cool  fluid  comes  into  contact  with 
atmospheric  air,  water  vapor  in  the  air  may  condense  or 
solidify.  This  is  the  problem  of  condensation/icing. 

Each  of  these  considerations  is  discussed  in  detail  in 
the  paragraphs  which  follow. 


Momentum  Mixing 


To  determine  the  effect  of  a  velocity  discontinuity 
between  the  cool  fluid  and  the  tunnel  fluid  it  is  con¬ 
venient  to  study  the  flow  situation  of  fig  2.  Two  fluids, 
both  of  which  have  infinite  extent  in  the  y-direction, 
meet  at  a  point  where  their  velocities  are  different.  The 
shear  stress  created  at  this  discontinuity  will  act  to 
smooth  the  discontinuity  and  a  momentum  mixing  region  will 
develop . 

This  flow  situation  is  discussed  in  detail  by 
Schlichting  (ref  1  )  and  by  Yih  (ref  5).  From  ref  1  , 
the  velocity  profile  in  the  mixing  region  is  given  by 

U ( y , t )  =  £ (Ui  +  U2)  +  £(Ui  -  U2)  X  {(3/2b)y  -  ( 1 /2b3 ) y 3 } 


The  quantity  b  which  appears  in  equation  17  is 
a  measure  of  the  size  of  the  momentum  mixing  region,  and 
Schlichting  gives 


b  =  3/2  6  2  (Ui  -  U2)t 


where  3  is  an  empirical  constant  which  depends  on  the 
(molecular)  mixing  length  and  must  be  determined  from 
experiment . 

The  important  conclusion  to  be  drawn  from  equation 
18  is  the  fact  that  the  mixing  region  will  grow  with  time 

The  flow  situation  for  which  high  Reynolds  number 
simulation  is  to  be  attempted  is  somewhat  different  from 
the  situation  depicted  in  fig  2.  This  is  because  the 
structure  over  which  the  cool  fluid  is  to  flow  limits  the 
extent  of  the  momentum  mixing  region. 

Equation  18  suggests  that  if  the  cool  fluid  velocity 
differs  from  the  tunnel  free  stream  velocity,  momentum 
mixing  will  eventually  force  the  two  velocities  to  be 
equal  in  the  case  where  the  cool  fluid  is  in  the  vicinity 
of  a  solid  boundary. 

This  equality  of  velocities  is  a  steady-state 
condition.  The  time  required  to  reach  steady-state  is 
not  easily  determined  analytically,  since  such  a  determ¬ 
ination  would  require  advance  knowledge  of  the  values  for 
8,  Ui,  and  U2.  Some  or  all  of  these  values  can  only  be 
determined  experimentally. 


wsv 


Nonetheless,  it  seems  reasonable  to  expect  that  the 
time  to  reach  steady-state  will  be  short.  This  is  because 
the  height  of  the  cool  fluid  layer  above  the  solid  boundary 
will  be  very  small  (on  the  order  of  magnitude  of  the  momen¬ 
tum  boundary  layer  height)  so  that  it  will  be  readily  influ 
enced  by  the  tunnel  free  stream. 

The  steady-state  assumption  should  be  confirmed  exper¬ 
imentally,  but  it  will  be  made  here  in  order  to  proceed  wit! 
the  analysis  of  the  proposed  high  Reynolds  number  sim¬ 
ulation  technique. 

As  a  result  of  the  discussion  above,  the  following 
result  is  indicated: 

"In  determining  the  parameters  required  for 
successful  high  Reynolds  number  simulation, 
the  free  stream  velocity  of  the  cool  fluid 
may  be  taken  as  equal  to  the  free  stream 
velocity  of  the  wind  tunnel." 


Thermal  Mixing 

The  success  of  the  high  Reynolds  number  simulation 
technique  under  consideration  depends  on  the  fact  that 
the  kinematic  viscosity  of  a  gas  decreases  with  decreasing 
temperature.  It  is  therefore  important  to  insure  that 
the  cool  fluid  stays  cool  and  is  not  significantly  heated 


by  the  much  warmer  tunnel  fluid. 

The  flow  situation  of  interest  is  depicted  in  fig  3. 


Appendix  B  to  this  report  contains  a  detailed  analysis 
of  this  situation  and  develops  the  solution  for  the  temp¬ 
erature  distribution  in  the  thermal  mixing  region. 

Since  the  energy  conducted  at  the  boundary  of  the 
thermal  mixing  region  can  be  balanced  by  the  convection 
of  fluid  into  the  thermal  mixing  region,  a  steady-state 
condition  will  occur  in  which  the  width  of  the  mixing 
region  will  not  grow  with  time. 

Under  the  assumptions  of  1)  steady-state;  2)  the 
temperature  mixing  region  does  not  effect  the  velocity 
field;  and  3)  gradients  in  the  y-direction  are  much  larger 
than  gradients  in  the  x-direction,  the  temperature  profile 
in  the  thermal  mixing  region  is  given  by 

0(x,y)  =  £erfc(y//4ax/V~)  (y  >  0)  (19) 

9 ( x,y)  =  1  -  serf c( | y | //4ax/ V^)  (y  <  0)  (20) 

where  0  is  a  non-dimensional  temperature  difference  defined 
by  0  =  (T  -  T  )/(T  -T). 

Of  the  assumptions  used  in  developing  these  expressions 
the  first  is  certainly  reasonable,  based  on  the  argument 
given  above.  The  third  assumption  is  just  a  restatement  of 
the  "boundary  layer  assumption,"  and  it  seems  almost  self- 
evident  that  the  thermal  mixing  region  will  have  a  boundary 
layer  characterisitic  to  it. 

The  second  assumption,  however,  is  not  quite  so  self- 
evident.  It  should  be  noted  that  Yih  (ref  5  )  makes  a 


similar  assumption  in  dealing  with  the  temperature  dis¬ 
tribution  created  by  a  pre-heated  jet. 

If  the  temperature  distribution  in  the  thermal  mixing 
region  were  to  disturb  the  velocity  field,  it  would  create 
a  situation  where  momentum  mixing  would  occur  in  the  free 
stream.  This  mixing  would  eventually  effect  the  develop¬ 
ment  of  the  momentum  boundary  layer  and  the  validity  of  the 
Reynolds  number  simulation  technique  would  be  question¬ 
able. 

Fortunately,  results  presented  below  suggest  that  the 
size  of  the  thermal  mixing  region  is  very  small  -  on  the 
order  of  magnitude  of  the  momentum  boundary  layer.  It 
seems  reasonable  to  assume  that  such  a  small  region  will 
not  have  a  significant  effect  on  the  velocity  field. 

For  the  sake  of  establishing  the  width  of  the  thermal 
mixing  region  it  is  necessary  to  determine  the  distance 
between  the  isotherms  0=1  and  0=0.  However,  as  is 
characteristic  in  boundary  layer  work,  the  solution  for 
0  indicates  that  the  thermal  mixing  region  extends  to 
infinity  in  both  directions.  Therefore,  we  adopt  the  con¬ 
vention  that  the  width  of  the  thermal  mixing  region,  2b^ , 
is  given  by  the  distance  between  the  isotherms  0  =  0.99 
and  0  =  0.01.  If  y  =  0  is  taken  as  the  line  of  symmetry 
between  these  isotherms,  then  b^  is  the  distance  from  y  =  0 
to  0  =  0.99. 

For  the  purpose  of  establishing  a  Reynolds  number 
simulation  technique,  the  extent  of  the  Lhermal  mixing 
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region  into  the  cool  fluid  must  be  studied.  In  terms  of 
the  flow  situation  depicted  in  fig  3  this  is  the  region 
y  <  0,  so  equation  20  yields 

9(x,y)  =  1  -  serf  c(  |  y  |  /  /Totx/V^)  =  0.99  (21) 

From  this,  erf c(  |y]//4ax/Voo)  =  0.02,  which  yields  the 
result  y//4ax/V  =  -1.65.  Since  b,  is  given  by  the  distance 
from  the  curve  y  =  0  to  the  curve  0  =  0.99,  it  is  clear  that 


Equation  22  does  not  apply  exactly  to  flow  over  a 
solid  boundary.  The  fluid  streamlines  will  be  deflected 
away  from  the  boundary  due  to  the  action  of  viscous  forces, 
so  the  isotherm  0  =  0.99  will  be  deflected  away  from  the 
boundary  as  well.  Therefore,  equation  22  can  be  taken  as 
a  conservative  estimate  for  the  extent  of  the  thermal  mixing 
region  into  the  cool  fluid  layer. 

If  the  thermal  mixing  region  extends  into  the  mom¬ 
entum  boundary  layer  it  will  effect  the  boundary  layer 
development.  A  condition  for  accurate  Reynolds  simulation 
is  therefore 

"Flow  properties  must  be  established  to 
guarantee  that  the  thermal  mixing  region 
does  not  extend  into  the  momentum  boundary 


Condensation  and  Icing 


From  the  discussion  of  the  preceding  paragraphs  it 
is  apparent  that  it  is  possible  to  create  a  flow  situation 
in  which  the  fluid  is  kept  separated  into  cool  and  warm 
regions.  If  the  warmer  fluid  contains  water  vapor  (as 
would  be  the  case  for  an  atmospheric  wind  tunnel)  the 
water  vapor  contained  within  the  thermal  mixing  region  can 
be  expected  to  cendense  and/or  solidify. 

This  presents  a  difficulty  in  attempting  to  simulate 
high  Reynolds  number  through  injection  of  a  cool  gas,  since 
the  water  and  ice  particles  will  be  significantly  more  dense 
than  the  surrounding  fluid.  The  motion  of  these  particles 
will  therefore  be  influenced  by  gravitational  acceleration; 
this  possibility  must  be  accounted  for. 

Water  and  ice  must  be  prevented  from  entering  the 
boundary  layer  if  accurate  Reynolds  number  simulation  is 
to  be  obtained.  Appendix  C  to  this  report  presents  a 
method  for  determining  the  trajectory  of  a  spherical 
particle  of  liquid  or  ice.  In  the  case  where  the  free 
stream  fluid  velocity  is  perpendicular  to  the  gravita¬ 
tional  field  the  equation  of  motion  for  the  sphere  is 
given  by 


v  =  -g*  +  C^(  3/ 8r)  pv2  if  |  g *  |  >  |  C<D(  3/8r )  pv2  | 
v  =  0,  otherwise  (2?) 
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subject  to  the  initial  condition  v(0)  =  0.  The  component 
velocities  of  the  particle  tangential  and  normal  to  the 
free  stream  are  denoted  by  u(t)  and  v(t),  respectively; 
p  is  the  ratio  of  the  free  stream  fluid  density  to  the 
density  of  the  particle  (p  =  Pf/PD)»  an<^  g*  as  the  effective 
acceleration  of  gravity,  taking  into  account  buoyant 
forces.  Cp  is  the  coefficient  of  drag  for  a  sphere,  taken 
to  be 

CD  =  24/Re  ( 0<Re<  0.25) 

5  =  0.102C2  -  0.894C  +  1-42  ( 0 . 25<Re<2000)  (24) 

CD  -  0.405  (2000<Re<200,000) 

where  z,  -  log^g  C^,  £  -  log^  Q  Re,  and  Re  is  based  on  the 
velocity  component  v. 

With  these  relationships,  equation  23  can  be  integ¬ 
rated  numerically  to  determine  particle  trajectory  as  a 
function  of  free  stream  velocity  and  sphere  diameter. 

Table  1  gives  the  distance  that  a  sphere  with  a  diameter 
of  approximately  V0  inch  can  be  expected  to  fall  as  a 
function  of  time.  Larger  particles  will  drop  more  quickly, 
smaller  ones  will  not  fall  as  fast. 

Since  the  mass  fraction  of  water  vapor  in  the  atmospher 
is  very  small,  any  condensation  which  occurs  in  the  tunnel 
can  be  expected  to  take  the  form  of  a  mist.  Therefore, 
large  particles  of  water  and  ice  (on  the  order  of  Vs  inch, 
for  instance)  will  not  form  in  the  tunnel.  Since  the 


results  presented  in  Table  1  and  in  Appendix  C  indicate 
that  even  relatively  large  particles  will  not  effect  the 
momentum  boundary  layer  for  reasonable  values  of  the  free 
stream  tunnel  velocity,  it  is  clear  that  any  water  vapor 
mist  will  be  insignificant,  as  well. 

Of  course,  it  is  also  possible  to  align  a  model  to 
be  tested  in  a  wind  tunnel  in  such  a  way  as  to  prevent  ice 
and  water  from  falling  towards  the  boundary  layer  (i.e., 
the  model  could  be  mounted  sideways) .  In  some  cases  for 
which  the  flow  geometry  is  not  as  simple  as  the  one  dis¬ 
cussed  above,  this  may  be  the  only  way  to  prevent  water 
and  ice  from  effecting  the  Reynolds  simulation  technique. 

The  conclusion  of  this  discussion  on  condensation  and 
icing  is  therefore 

"  Condensation  and  icing  must  be  prevented 
from  effecting  the  momentum  boundary  layer 
development  either  by  guaranteeing  that 
fluid  free  stream  velocity  is  large  enough 
to  sweep  moisture  away  from  the  region  of 
concern,  or  by  choosing  an  appropriate 
orientation  for  the  model." 

Other  Considerations 

It  is  clear  that  a  primary  consideration  in  the 
development  of  a  successful  Reynolds  number  simulation 
technique  is  the  necessity  to  prevent  mixing  between  the 
warm  and  the  cool  fluids .  The  discussion  above  outlines 
ways  of  dealing  with  this  concern  for  the  flow  phenomena 


which  will  effect  mixing  most  significantly. 


Three  additional  considerations  may  impact  mixing, 
as  well.  If  the  cool  fluid  is  not  the  same  as  the  warm 
fluid,  diffusion  will  occur  between  the  two  fluids.  There¬ 
fore,  it  is  necessary  to  select  fluids  which  are  "similar 
enough"  (air  and  nitrogen,  for  example)  to  minimize  the 
effect  of  this  diffusion. 

The  static  pressure  of  the  two  fluids  must  be  the 
same  in  the  region  of  interest.  This  is  not  expected  to 
be  a  significant  problem  for  the  technique  under  consid¬ 
eration,  since  the  cool  fluid  will  be  injected  through 
some  type  of  nozzle.  When  the  flow  through  the  nozzle  is 
subsonic,  the  fluid  leaving  the  nozzle  will  have  a  static 
pressure  equal  to  the  static  pressure  of  the  tunnel  fluid. 

The  last  consideration  is  the  development  of  density 
gradients.  Gravity  will  pull  a  more  dense  fluid  down 
into  a  less  dense  fluid,  and,  in  general,  the  cool  fluid 
will  be  more  dense  than  the  boundary  layer.  In  the  case 
of  an  attempt  to  simulate  boundary  layer  development  on 
the  top  surface  of  some  object,  the  density  gradient 
should  not  be  a  significant  concern  since  the  cool  fluid 
will  be  below  the  warm  fluid.  Across  the  bottom  surface, 
the  density  gradient  will  be  significant  unless  the  inertia 
of  the  fluid  is  sufficient  to  overwhelm  the  effect  of  the 
gravitational  force  (i.e.,  when  Gr/Pe2  <<  l).  In  most 


other  cases  the  model  to  be  examined  will  have  to  be 


h  Reynolds  Number  Simulation:  the  Flat  Plate 


This  section  considers  in  detail  Reynolds  number 
simulation  through  injection  of  a  cool  gas  for  the  case 
of  a  flat  plate. 

The  flow  situation  of  interest  is  illustrated  in  fig 
4.  A  fluid  with  uniform  velocity,  but  separated  into 
layers  of  cool  and  warm  temperature,  approaches  the  leading 
edge  of  a  flat  plate  located  at  position  x  =  0.  The 
x-axis  is  aligned  with  the  interface  between  the  two  tem¬ 
perature  regions  and  is  located  at  a  height  h  above  the 
flat  plate.  The  momentum  boundary  layer  thickness  at  a 
point  x  is  denoted  by  6 ,  and  2b^  gives  the  width  of  the 
thermal  mixing  region.  y^  and  y^  are  respectively  the 
y-coordinate  of  the  lower  boundary  of  the  thermal  mixing 
region  (the  curve  0  =  0.99)  and  the  ordinate  of  the  momen¬ 
tum  boundary  layer  thickness.  Both  of  y^  and  y^  are  func¬ 
tions  of  x. 

The  point  at  which  y^  =  y^  has  been  denoted  by 
x  ,  it  is  at  this  location  that  the  thermal  mixing 
region  will  begin  to  effect  the  development  of  the  momen¬ 
tum  boundary  layer. 

T^  is  the  temperature  of  the  flat  plate,  and  for 
the  purpose  of  this  study  it  is  assumed  that  uo  heat  trans¬ 
fer  is  taking  place  at  the  plate.  Therefore,  Tw  is 
assumed  equal  to  Tc  (or  to  the  adiabatic  wall  temperature 


for  the  case  of  high  speed  flow) . 
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It  has  been  established  in  previous  chapters  that 
a  necessary  condition  for  successful  Reynolds  number  sim¬ 
ulation  prohibits  the  thermal  mixing  region  from  entering 
the  momentum  boundary  layer. 

For  a  flat  plate  of  length  t,  this  condition  requires 

that  x  . ,  >  t. 
cnt 

For  laminar  flow,  Schlichting  (ref  1  )  gives  the 
boundary  layer  height,  6,  as 

6  -  5.0/vx/V^  (25 

In  the  coordinate  reference  frame  defined  by  fig  4»  this 

results  in  the  relationship  y^  =  -h  +  6.  It  has  already 

been  established  that  =  -3 . 3/otx/V^ .  The  abscissa  for 

which  y,  =  y  is  x  so 

Jt  Jm  crit' 

-3.3/o Tpr /x  77  =  -h  +  5 . 0/v/ V  /x  77  (26 

°°  crit  «>  cnt 

where  the  values  of  a  and  v  are  determined  at  T  . 

c 

Solving  equation  26  for  x  yields  the  result 

xcrit  =  127 

where  is  defined  by  L ^  =  5 . 0/v/ +  3 . 3/a/ Vro .  The 

condition  x  . ,  >  L  then  forces  the  result 
crit 


(k/L,)2  >  l 


(28 


\  V 


which,  as  a  constraint  on  h,  is  written  as 


h  >  L^/Z 


A  similar  condition  can  be  derived  for  the  case  of 
turbulent  flow  by  defining  the  parameter  L For  turbulent 
flow  over  a  flat  plate,  Holman  (ref  6)  gives 

6  =  0.37 x^5  (30) 

This  value  for  6  is  substituted  into  the  expression  ym  = 

-h  +  6  and  y^  is  equated  to  ,  as  above.  As  before,  this 
yields  the  result 


xcrit  =  WLt) 


where  L ^  is  defined  by  L  ^  -  3 . 3/ot/ Vro  +  0 . 37^v/Vjxc1°  ^  . 
Applying  the  constraint  xcr^  >  l  and  defining  1 ^  in 
terms  of  the  length  £  (rather  than  xcr:j^)  yields  the 


condition 


h  >  L  jJZ 


Equations  29  and  32  are  the  mathematical  expressions 
of  the  requirement  that  the  thermal  mixing  region  may  not 
extend  into  the  momentum  boundary  layer.  There  are,  how¬ 
ever,  additional  constraints  on  the  value  of  h. 


One  apparent  constraint  on  h  is  the  requirement  that 


h  >  5,  to  guarantee  that  the  momentum  boundary  layer  is 
contained  within  the  cool  fluid.  This  condition  is  auto¬ 
matically  satisfied  if  k  >  L/Z,  since  L/Z  >  6  in  all  cases. 

Another  constraint  on  h  is  related  to  the  mass  rate  of 
flow  of  the  cool  fluid,  m.  The  mass  rate  of  flow  is 
directly  proportional  to  the  cross-sectional  area  across 
which  the  fluid  moves.  For  the  case  of  the  flat  plate, 
this  cross  sectional  area  is  given  by  hx u),  where  w  is 
the  width  of  the  plate.  Therefore, 


-  =  p  V  h 
M  M  c  00 


(33) 


where  pc  is  the  density  of  the  cool  fluid. 

Since  it  has  been  established  above  that  h  >  L/Z, 
equation  33  can  be  cast  as  a  constraint  on  m. 

—  >  p  V  L/Z  (34) 

Therefore,  given  the  value  of  it  is  possible  to 

determine  the  required  minimum  value  of  m. 

If  the  cool  fluid  can  be  treated  as  an  ideal  gas, 

then  p  can  be  determined  from  T  and  p  .  It  has 
c  c  ^c 

already  been  shown  that  pc  must  equal  the  static  pressure 

(p  )  of  the  warm  fluid.  Therefore,  it  is  assumed  that  p 
s  c 

p  ,  and  then  p  can  be  determined  from  the  properties  of 
s  s 


the  warm  fluid. 


If  the  flow  of  the  warm  fluid  can  be  assumed  to  be 
isentropic,  the  following  algorithm  can  be  used  to  determine 


1)  Given  V  ,  T,  ,  and  p, ,  guess  p  . 

00  X*  u  S 

2)  Determine  q  =  IpV^2  . 

3)  Let  ps  =  pt  -  q 

4)  Form  the  ratio  (p  /p ,  )  ,  and  apply  the 

S  b 

isentropic  relationships  to  determine 
the  ratio  (T  /T, ): 

S  w 

a)  Guess  a  value  for  y. 

b)  Solve  Tg  =  Tt(ps/pt)(Y-1)/^ 

c)  Determine  y'  =  y'(T  ) 

s 

d)  If  y'  =  Y  then  stop,  otherwise  y  =  y'. 
Repeat  steps  b  through  d. 

5)  Determine  p  *  =  p  /(RT  ). 

s  s  s 

6)  If  p  '  =  p  ,  then  stop.  Otherwise, 

s  s 

let  p  =  p  '  and  iterate  steps  2 

b  b 

through  6. 


In  many  cases,  the  assumption  that  p  =  p,  will 
be  reasonably  accurate  (see  Appendix  D) ,  so  it  may  not 
be  necessary  to  employ  the  algorithm  defined  above  in 
order  to  obtain  quick  estimates  for  fluid  properties. 

Once  p  has  been  determined,  the  value  of  p  can 
be  found  from  the  ideal  gas  law.  Since  V  is  assumed 
equal  to  the  free  stream  tunnel  velocity,  it  is  possible 
to  calculate  the  required  mass  rate  of  flow  of  the  cool 
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fluid  from  equation  34. 

A  final  constraint  on  the  value  of  h  involves  the 
concern  over  possible  condensation  and  icing.  In  the  worst 
case  (assuming,  again,  that  there  is  no  momentum  mixing 
occurring)  ice  or  water  may  form  beginning  at  the  point 
where  the  two  layers  of  fluid  are  allowed  to  exchange  energy 
In  the  notation  given  in  fig  4  this  is  the  point  (x,y)  =  (0, 

A  spherical  particle  formed  at  location  (0,0)  and 
moving  with  velocity  in  the  x-direction  will  traverse 
the  length  of  the  plate  in  a  period  of  time  £/V  .  The 
equations  of  motion  established  in  Appendix  C  to  this  report 
can  be  used  to  determine  the  distance  that  a  particle  of 
given  diameter  will  fall  in  the  time  period  The 

height  of  the  cool  fluid  layer  can  then  be  selected  to 
insure  that  the  particle  will  not  fall  into  the  boundary 
layer . 

In  practice,  it  is  probably  easiest  to  select  h  and 
m  based  on  the  criteria  given  in  equations  29  or  32  and 
34  and  then  to  check  and  see  if  this  height  is  sufficient 
to  prevent  water  and  ice  from  falling  into  the  boundary 
layer.  If  not,  h  and  m  may  be  adjusted  accordingly. 

The  results  presented  above  may  be  simplified  some¬ 
what  if  cruder  estimates  for  the  values  of  the  important 
parameters  are  acceptable. 

From  the  definition  of  ,  it  is  apparent  that 
~  26  f  since  a  ~  v  for  most  gases.  Therefore,  the  con- 
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straint  on  h  can  be  expressed  approximately  a: 


h  >  26  (35) 


and  the  constraint  on  the  mass  rate  of  flow  becomes 


-  >  2p  V  5  (36) 

M  *C  oo  ' 


Ice  or  water  should  not  effect  the  boundary  layer 
development  if  the  condensation  does  not  enter  the  boun¬ 
dary  layer.  Since  h  ~  26,  this  requires  that  the  distance 

a  solid  or  liquid  particle  may  fall,  y  ,  should  be  less 

P 

than  h/2.  If  viscous  forces  on  the  particle  are  neglec¬ 
ted,  its  equation  of  motion  can  be  approximated  by 


y p  -  sgt2 


The  time  t  is  given  by  -6/V^,  and  g  32  ft/sec2 
In  this  case,  equation  37  becomes 


(37) 


v  </ 

Vi*.- 


.  -  .  " 

V  V 


y  p  *  i6(^/vj2 


The  constraint  on  y^  is  given  by  y^  <  h/2 ,  so 


(38) 


16 (l/V^)2  <  or 


32 (£/Vj2<  lt/Z 


(39) 


which  can  be  written  in  terms  of  6  as 


16 U/V  )2  <  6 


For  turbulent  flows,  the  momentum  boundary  layer 
thickness  is  generally  much  greater  than  the  depth  of  the 
thermal  mixing  region,  so  ~  6.  In  this  case,  the 

constraints  on  h  and  m  can  be  approximated  by 


h  >  6 


->  p  V  6 

w  M  c  00 


Direct  application  of  equations  41  and  42  would 
result  in  the  constraint  y^  <  0,  which  is  clearly  unaccept 
able.  It  would  be  convenient  to  develop  a  relationship 
similar  to  equation  39  for  the  case  of  turbulent  flows. 

It  can  be  argued  that  the  high  energy  characteristic  of 
turbulent  flows  will  tend  to  overwhelm  the  effects  of 
condensation  in  the  boundary  layer,  so  some  penetration 
of  the  momentum  boundary  layer  by  water  or  ice  may  be 
tolerated.  Guided  by  the  form  of  equation  39,  the  con¬ 
straint  on  y^  is  again  chosen  to  be  y^  <  h/2  to  yield 


32 U/VJ2  <  Lt/Z 


In  summary,  the  constraints  which  determine  the 
necessary  parametric  values  for  successful  Reynolds  number 


simulation  are: 


1  .  h  >  L/Z 


p  ifi 

'  w 


P  V 

C  c 


Li/Z 


3.  32 U/VJ2  <  Lv>T 


Some  results  regarding  the  parameter  L  which  appear 
in  each  of  these  expressions  are  developed  in  Appendix  7 
to  this  report. 


VI.  A  Proposed  Experimental  Design 

For  the  purpose  of  verifying  some  aspects  of  the 
theoretical  development  presented  in  this  study,  the 
experimental  apparatus  depicted  in  fig  5  has  been  designed 
and  built  with  the  assistance  of  the  AFIT  shop. 

The  apparatus  is  basically  a  thick  flat  plate  made 
of  aluminum.  An  ellipsoidal  nose  of  plexiglas  is  designed 
to  fit  on  the  leading  edge  of  the  plate,  the  cool  fluid 
is  to  be  injected  through  this  nose  and  across  the  plate. 

The  plate  is  one  inch  thick  and  is  constructed  from 
two  one-half  inch  thick  sheets  of  aluminum  which  are 
bolted  together.  Three  small  channels  lead  from  a  mani¬ 
fold  located  near  the  leading  edge  of  the  plate  to  open¬ 
ings  at  the  trailing  edge  so  that  coolant  may  be  fed  to 
the  inside  of  the  plate. 

The  necessity  of  these  cooling  channels  is  demon¬ 
strated  by  a  lumped-capacity  analysis  of  the  plate.  In 
a  lumped-capacity  analysis,  the  temperature  of  the  alum¬ 
inum  is  assumed  to  be  uniform  and  its  dependence  on  time 
is  given  by 

T  -  T 

m - Tp-  =  exp{  -  (hs/pc  V)t)  (44) 

i  0  -  1  P 

u  no  1 


where  T0  is  the  initial  temperature  of  the  plate  and 
T  is  the  temnerature  of  the  free  stream  fluid  around  the 

OO  A 


plate,  V  is  the  volume  of  the  plate  and  S  is  its  surface 
area,  h  is  the  convective  heat  transfer  coefficient,  and 
p  and  0^  are  the  density  and  specific  heat,  respectively 
of  the  plate. 

Equation  44  is  a  good  approximation  for  T(t)  when 


the  Biot  number,  Bi 


■,  is  much  less  than  one, 


flat  plate  with  thickness  a,  the  Biot  number  is  given  by 


Bi  = 


The  thermal  conductivity  of  aluminum  can  be  taken 


to  be  240  W/m-K  and  the  value  of  the  heat  transfer  coef¬ 
ficient  can  be  obtained  from  the  empirical  relationship 

=  Nu  =  0.0296  Re^5  Pr^3  for  turbu  Lent  flow  (ref  6). 
a 

For  the  plate  under  consideration,  1-2  ft,  and 

k  ,  the  thermal  conductivity  of  air,  can  be  taken  to  be 
a 

0.03  W/m-K. 


With  these  values,  and  for  h  evaluated  at  x  =  l , 
the  value  of  Bi  can  be  obtained  as  a  function  of  free 
stream  fluid  velocity.  In  the  present  case,  0.0015  <  BI 
<  0.0168  for  15  m/sec  <  <  305  m/sec,  so  a  lumped  cap¬ 

acity  analysis  is  justified. 

Denoting  by  t  the  quantity  (pc.V/hS)  found  in 
c  p 


equation  44,  the  time  to  reach  steady- stat-. 


;mr  a  t,  u  r  e 


can  be  chosen  (by  convention)  to  be  3 1  .  For  the  veloc¬ 
ity  range  noted  above,  this  yields  660  sec  >  '  _  s  eU  s e 0 
These  values  indicate  an  unacceptably  lor.  •  per •  od  f  1  hr 
for  the  establishment  of  steady-state  coni',*  '  rs.  r 

this  reason  (and  because  it  has  been  assures  1  ha  *  '  : 


is  no  heat  transfer  at  the  plate)  it  is  necessary  to 
provide  a  means  of  cooling  the  plate. 

The  width  of  the  plate  is  one  foot  so  that  ’  1  may  be 
tested  in  the  AFIT  14-inch  wind  tunnel.  Ten  static  pres¬ 
sure  taps  are  located  4-5  in  from  one  side  of  the  plate  and 
are  evenly  spaced  down  its  length.  Five  copper-constantan 
thermocouples  are  situated  4.5  in  from  the  other  side  of 
the  plate. 

The  plexiglas  nose  is  designed  to  have  a  uniformly 
varying  radius  of  curvature  so  that  the  fluid  moving  over 
it  will  not  experience  an  abrupt  change  in  radius  of 
curvature  as  it  reaches  the  leading  edge  of  the  plate.  The 
equation  for  the  nose  is  given  by  (— )n  +  (r)n  =  1  with 
a  =  1.25  in,  b  =  0.5  in,  and  n  =  2.5.  The  nose  is  13  in  wide 
and  has  a  0.5  in  hole  drilled  through  it.  Gold  nitrogen 
can  be  fed  to  this  channel  as  either  a  liquid  or  as  a  gas. 

A  slot  of  0.03  in  is  cut  at  the  top  of  the  channel 
at  the  shoulder  point  of  the  ellipsoid  (see  fig  5  ).  The 
top  of  the  nose  is  cut  smoothly  to  meet  the  slot,  and  is 
designed  to  help  turn  the  nitrogen  so  that  it  will  flow 
across  the  plate  and  parallel  to  the  tunnel  free  stream. 

The  nose  is  designed  with  a  slot  that  fits  a  groove 
on  the  leading  edge  of  the  aluminum  plate.  This  allows 
the  possibility  of  testing  a  variety  of  different  nozzle 
designs  by  simply  interchanging  plexiglas  noses. 

Cryogenic  nitrogen  is  chosen  for  use  as  the  cool 
fluid  due  to  its  similarity  to  atmospheric  air  (i.c.. 


atmospheric  air  is  approximately  78%  nitrogen).  As 

noted  before,  the  nitrogen  can  be  introduced  into  the 

nose  as  either  a  liquid  or  a  gas;  the  choice  of  phase  is 

dependent  on  the  desired  values  for  T  and  m. 
r  c 

Use  of  gaseous  nitrogen  allows  control  over  the 
pressure  inside  the  nozzle,  and  therefore  ft  can  be  easily 
controlled.  A  disadvantage  to  the  use  of  gaseous  nit¬ 
rogen,  however,  is  the  difficulty  in  controlling  and  mon¬ 
itoring  its  temperature. 

The  use  of  liquid  nitrogen  is  constrained  by  a  limited 
ability  to  control  m  -  the  mass  rate  of  flow  will  be  prim¬ 
arily  dependent  on  the  rate  of  evaporation  of  the  nitrogen 
within  the  nozzle.  This  rate  can  be  controlled  to  a  lim¬ 
ited  degree  by  providing  a  means  of  heating  the  plexiglas 
nose . 

Liquid  nitrogen  could  also  be  boiled  externally,  and 
then  fed  to  the  nose  as  a  gas.  This  external  boiling 
would  all^w  greater  control  over  ft,  and  T  could  be  con- 
trolled  by  providing  sufficient  insulation  for  the  fluid. 

For  the  purpose  of  studying  the  theoretical  pred¬ 
ictions  of  this  paper  there  is  an  additional  advantage  to 
the  use  of  liquid  nitrogen.  Since  the  temperature  inside 
the  nozzle  will  be  the  boiling  temperature  of  nitrogen 
(~  140  R)  and  since  the  expansion  through  the  nozzle  will 
be  isenthalpic  (a  Joule-Thomson  expansion),  the  temperature 
of  the  gas  leaving  the  nozzle  will  be  known. 

Measurement  of  the  flow  properties  presents  a  chal- 


lenge  because  of  the  very  low  temperatures  involved. 

Dynamic  pressure  and  temperature  must  be  measured  simul¬ 
taneously  in  order  to  determine  the  velocity  profile  in 
the  boundary  layer.  The  dynamic  pressure,  q,  is  equal  to 
2 p v 2 ,  and  the  temperature  is  required  to  evaluate  p.  Tem¬ 
perature  readings  are  also  required  in  order  to  compare 
the  predicted  extent  of  the  thermal  mixing  region  with 
uieasured  data. 

For  accurate  measurements  with  reasonable  response 

time,  hot  wire  or  hot  film  anemometers  may  be  used.  It 

should  be  noted  that  hot  wire  (film)  probes  are  calibrated 

for  a  particular  temperature;  if  they  are  not  calibrated 

at  T  a  correction  must  be  made  for  this.  This  is  dis- 
c 

cussed  in  some  detail  in  ref  7,  and  ref  8  gives  a  method 
for  using  hot  wire  anemometers  to  measure  both  temperature 
and  dynamic  pressure  under  very  controlled  conditions. 

For  a  quick  check  on  the  predicitions  made  in  this 
paper  it  may  suffice  to  measure  q  with  a  pitot  probe  and 
to  take  total  temperature  readings  with  a  thermocouple. 

If  the  flow  speed  is  not  too  high,  the  total  temperature 
will  be  an  adequate  approximation  to  the  static  temperature 
of  the  fluid. 

Care  must  be  taken  at  the  beginning  of  the  test  run 
to  insure  that  the  surface  of  the  plate  does  not  ice  - 
the  nitrogen  must  be  swept  quickly  across  the  plate  to 
displace  the  air  before  water  vapor  in  the  air  can  con¬ 


dense  . 


This  difficulty  can  be  eliminated  altogether  if  warm 
nitrogen  is  used  to  sweep  the  atmospheric  air  from  the 
surface  of  the  plate  before  the  cool  nitrogen  is  introduced 

The  analyses  contained  in  Appendices  D  and  E  to  this 
report  indicate  that  very  high  mass  rates  of  flow  for  cool 
nitrogen  are  required  for  high  Reynolds  number  simulation. 
The  experimental  apparatus  described  above  will  not  be  able 
to  handle  those  requirements,  due  primarily  to  the  small 
size  of  the  nozzle  through  which  the  nitrogen  is  being 
injected. 

Nonetheless,  this  apparatus  can  be  used  to  provide 
preliminary  verification  of  certain  aspects  of  the  theor¬ 
etical  development  given  in  this  paper.  In  particular, 
experimental  data  is  required  to  more  accurately  determine 
the  extent  of  the  thermal  mixing  region,  the  effect  of  any 
momentum  mixing  which  may  occur,  and  the  impact  (if  any) 
of  condensation  or  icing  on  momentum  boundary  layer 
development . 


Conclusions  and  Recommendations 


This  paper  has  analyzed  the  potential  for  high  Rey¬ 
nolds  number  simulation  in  a  wind  tunnel  through  use  of  a 
cool  gas  to  envelop  the  momentum  boundary  layer.  If  this 
technique  can  be  used,  effective  Reynolds  number  simu¬ 
lation  can  theoretically  be  increased  by  a  factor  of  seven 
or  eight. 

To  be  successful,  the  flow  of  the  warm  tunnel  air 
must  not  interfere  with  the  locally  cool  fluid  flow.  For 
the  case  of  a  flat  plate,  it  has  been  shown  that  this 
condition  will  be  met  when: 

1  .  h  >  L/Z 

2.  m/w  >  p  V  L/Z 

c  00 

3.  32 (Z/yJ2  <  L/Z 

The  conditions  defined  above  are  based  in  part  on 
certain  assumptions  -  in  particular,  that  the  cool  fluid 
free  stream  velocity  is  equal  to  the  tunnel  free  stream 
velocity,  and  that  condensation  and  icing  are  not  sig¬ 
nificant  concerns  if  the  condensation  does  not  penetrate 
too  far  into  the  boundary  layer. 

Under  these  circumstances,  high  Reynolds  number  sim¬ 
ulation  can  theoretically  be  achieved  if  the  constraints 
1  -  3  noted  above  are  satisfied. 

In  practice,  the  required  value  of  ffi  for  high  Roy- 


\  \  ^  *.  N.  \  \  V'  *.  V  t  *.  \  '. 


nolds  nuraer  simulation  will  be  very  large  (on 
0.5  to  1.0  lb  /sec  for  the  example  discussed 
Therefore,  it  may  not  be  a  simple  matter  to  a 
etical  constraints  to  a  situation  of  interest 


%  \  A  •  -  -W  - 


the  order  of 
in  this  paper) 
pply  the  theor 


Recommendations  for  Further  Study 


Further  study  in  the  following  areas  will  add  signif¬ 
icantly  to  the  development  of  a  successful  high  Reynolds 
number  simulation  technique: 

1 .  The  experimental  apparatus  discussed  in  Chapter 
VI  should  be  used  to  verify  the  assumptions  noted  above. 

2.  The  development  which  has  been  presented  for  the 
case  of  a  flat  plate  should  be  extended  to  structures  of 
arbitrary  shape. 


3.  Various  means  of  injecting  the 
order  to  provide  the  mass  rates  of  flow 
successful  high  Reynolds  number  simulati 
investigated . 


cool  fluid  in 
required  for 
on  should  be 
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Appendix  A.  The  Governing  Equations  for  Viscous  Fluid  Flow 

The  purpose  of  this  appendix  is  to  develop  in  detail 
the  equations  presented  in  Chapter  II  which  model  the  behav¬ 
ior  of  viscous  fluids. 


The  Continuity  Equation 

The  mass  contained  in  a  control  volume,  ft,  at  any 
time,  t,  is  equal  to  the  volume  of  the  region  ft  multiplied 
by  the  average  density,  p(t),  of  the  mass  contained  within 
ft  at  that  time.  In  general,  mass  may  be  transported  across 
the  boundary  of  the  control  volume,  9ft.  The  excess  of 
the  mass  entering  the  region  over  the  mass  leaving  the  region 
must  remain  stored  within  ft. 

The  rate  of  mass  transfer  at  some  point  on  9ft  is 
given  by  (pV  )*  x  (As),  where  As  is  an  element  of  surface 
area  on  3ft,  Vn  is  the  component  of  the  fluid  velocity  per¬ 
pendicular  to  the  element  As,  and  ( p  V  r  )  ■'  is  the  average 

value  of  pV  taken  over  As.  For  convention,  pV  is  taken 
n  n 

as  positive  going  into  ft,  and  negative  otherwise. 

If  no  mass  is  created  within  ft,  the  mass  contained 
within  the  region  at  time  t+At  must  be  equal  to  the  mass 
stored  at  time  t  plus  the  net  influx  of  mass  into  the 
region  over  the  time  At.  Mathematically,  this  can  be 
expressed  as 


52 


p ( t+ At )  X  Vol(n)  =  p(t)  x  Vol(fi)  +  (pV  )*  X  Area(3Q)  x  At 

(A.1) 

Dividing  through  by  At,  rearranging,  and  taking  the  limit  as 
At  -x  0,  this  becomes 

||  x  Vol(fi)  =  (pV  )*  x  Area (3fi)  (A. 2) 

o  Xj  n 

As  Vol(Sl)  -x  dVol  so  Area(3ft)  ->•  ds,  and  average  properties 
taken  over  the  volume  Q  approach  the  average  values  taken 
over  3fi.  Furthermore,  pV^  can  be  rewritten  as  -(pv*n),  where 
n  is  taken  to  be  the  outward  pointing  normal  associated  with 
ds.  A  fundamental  result  from  the  calculus  of  vectors  and 
tensors  yields  pv*n  ds  =  V*pv  dVol  (ref  9  ).  Applying  these 
results  to  equation  A. 2  gives 

=  _  ( V«pv)  ( A . 3) 

This  is  the  Continuity  Equation,  and  it  is  more 
frequently  written  as 

+  (V*pv)  =  0,  or  equivalently  as 
+  p(V*v)  =0  (A.  4) 

In  many  applications,  it  is  reasonable  to  take  p 
as  essentially  constant,  in  which  case  ^  is  identically 
zero  and  the  Continuity  Equation  reduces  to 


The  Momentum  Equation 


The  momentum  contained  within  a  fixed  region  in  space 
( Q.)  will  change  with  time  due  to  the  addition  of  mass  to  fi 
and  due  to  forces  which  act  on  the  mass  contained  in  the 
control  volume.  Considering  the  limiting  case  as  Vol(ft)  -*  0 
this  can  be  expressed  mathematically  as 

(pv)(t+At)  dVol  =  pv(t)  dVol  +  (pvV^)  ds  At  +  (a  ds)  At 

+  (pf  dVol)  At  (A. 6) 

where  f  denotes  the  body  force  per  unit  mass  acting  on  the 
fluid  contained  within  and  o  is  the  force  due  to  stress 
acting  on  the  surface  3fl.  It  is  clear  that  o  depends  on 
the  state  of  stress  at  ds  as  well  as  on  the  orientation  of 
the  elemental  surface  area.  This  dependence  is  expressed 
as  o  =  o*n,  where  o  is  the  stress  tensor  (ref  10  )  . 

Noting  that  pvV^  =  -  pvvh,  dividing  through  the 
equation  above  by  At,  and  taking  the  limit  yields 

dVol  -  -  ^-pvv«nds  +  o*nds  +  pfdVol  (A.  7) 

At  equilibrium,  application  of  the  law  of  conservat,  i  on 
of  angular  momentum  demonstrates  that  o  is  symmetric.  fine* 


v  v  is  the  outer  product  of  a  vector  with  itself,  v  v 
also  symmetric,  and  so  o*n  ds  =  V*o  dVol  and  pv  v  • h  d 
V • ( pv  v  )  dVol.  The  term  V* ( pv  v  )  can  be  expanded  as 
(Vpv)v  +  p(vVv).1  The  continuity  equation  gives  Vi 
-  so  the  equation  above  reduces  to 


3  ( pv) 
3t 


-  v(-|^)  +  p(v*Vv)  =  Vo  +  pf 


Expanding  the  first  term  on  the  left  side  of  this  equation 
and  simplifying  yields 


+ 


v  •  Vv ) 

Dv 

pDt 


Vo  +  pf,  which  is  equivalent  to 
V*  a  +  pf 


(A. 


Equation  A. 9  is  Cauchy's  equation  of  motion.  The 
Navier-Stokes  equations  are  developed  from  equation  A. 9 
by  determining  the  relationship  between  the  stress  tensor 
and  the  rate  of  deformation  of  the  fluid  contained  within 
ft.  It  is  convenient  to  separate  the  stress  tensor  into 
components  perpendicular  to  3ft  (pressure)  and  tangential 
to  3ft  (shear  stress).  Since  pressure  is  taken  ns  constant 
over  3(9  for  ft  sufficiently  small,  wo  can  write  o  -  -pT  +  t 
where  I  is  the  identity  tensor. 

If  the  fluid  is  assumed  to  be  Newton :  nti ,  the  component 

1  Note  that  v*7v  i  (v*V)v  in  general  .  unless  a  a r4  •  s i an 
coord  i  nate  system  is  use!  .  A  general  i  dent.  Lty 
v  •  7  v  -  h  g  r  a  d  (  v  •  v )  -  (  v  x  (  V  *  v  )  ) 


%  \ 


of  stress  are  linearly  dependent  on  the  rates  of  deformatl 
of  the  fluid.  It  can  be  shown  (refs  1  and  11  )  that  this 
relationship  is  given  by  t  =  A(V*v)l  +  u(Vv  +  grad  v). 

Substituting  this  expression  for  t  into  equation  A. 9 
and  making  use  of  the  identity  V*al  =  Va  (a  is  a  scalar) 
results  in  the  Navier-St okes  equations: 

pjpp  =  -Vp  +  V(aV*v)  +  V*u(Vv  +  grad  v)  +  pf  (A. 

As  a  consequence  of  Stokes'  Hypothesis  (ref  12),  the 
parameters  u  and  A  are  assumed  related  by  the  equation 
3  A  +  2p  =  0. 

If  the  coefficient  of  viscosity,  p,  is  assumed  con¬ 
stant,  and  if  the  function  v  is  twice  continuously  differ¬ 
entiable  so  that  the  order  of  differentiation  may  be  inter 
changed  (resulting  in  the  identity  V*grad  v  =  V(V*v)  ), 
then  equation  A. 10  reduces  to 

p-jpp  -  -Vp  +  uV2v+V3u9(V*v)+pf  (A. 

In  the  case  of  incompressible  flow,  the  font  inu :  f y 
Equation  gives  V*v  =  0,  and  when  body  forces  a r •  -  1  '  - 

1  V2  v  =  V*Vv  /  (V«V)v,  except,  in  far4  -v  i  a . '•  •  ' 

s  y  s  t  e  m  s  . 

Alternately,  V2v  =  V(V*v)  -  (7-(,'-v  ‘ 


m 


p? 


r* 

N 


pi 


■v-r  ■■.  '•■ 


ted,  the  Momentum  Equation  takes  the  form 


Dv 


Dt 


=  -Vp  +  yV2v 


( A  .  1  2 ) 


The  Energy  Equation 


The  energy  contained  within  a  fixed  region  in  space 
(fi)  will  change  with  time  due  to  addition  of  mass  to  the 
control  volume,  heat  added  to  the  control  volume,  and  work 
done  on  the  control  volume.  Considering  the  limiting  case 
as  Vol(Sl)  -*■  0,  this  is  expressed  as 


pE(t+At)  dVol  =  pE(t)  dVol  +  (o*v)At  ds  +  pf*vAt  d; 


+  q  At  ds  +  (pvE)  At  d: 
Hn  n 


(A. 13) 


where  E  is  the  energy  per  unit  mass  and  q^  is  the  component 


of  heat  flux  perpendicular  to  ds  and  directed  into  fi.  By 


noting  that  o  =  o*n,  q^  =  -q*n,  and  (pvE)  =  -pvE*n,  this 


becomes 


pE(t+At)  dVol  =  pE(t)  dVol  +  (o*v)*nAt  ds  +  pf*vAt  ds 

-  q  •  n  At  ds  -  pvF,*nAt  ds  (A.  14) 


Applying  the  identity  a*h  ds  =  V*a  dVol  for  any 
vector  a,  dividing  through  by  At,  and  taking  the  limit  a; 


1  '/■ 


At  ■>  0,  the  energy  equation  takes  the  form 


— (pE)  =  V*(a*v)  +  pf»v  -  V*q  -  V*pvE 


(A.  1 


Applying  the  vector  identity  V*  ca  =  c(V»a)  +  a*Vc 
to  the  last  term  of  this  equation,  expanding  the  derivative 
on  the  left,  and  noting  that  =  -(V'p'O  from  continuity, 
the  energy  equation  takes  the  form 


V* ( a*  v )  +  pf • v 


V  •  q 


( A .  1 


This  is  the  total  energy  equation.  In  the  case 
where  changes  in  potential  and  chemical  energy  can  be 
neglected,  E  can  be  taken  to  be  (e  +  |v*v)  where  e  is 
the  internal  energy  per  unit  mass  and  gw  is  the  kinetic 
energy  per  unit  mass. 

When  the  momentum  equation  ( eqn  A. 9)  is  dotted  with 
the  velocity  vector,  the  result  is 

pjj^r«v  =  V»o*v  +  pf*v  (A.1 

When  this  is  subtracted  from  the  total  energy  equation 
(making  use  of  the  identity  V*(o*v)  =  V*o*v  +  o:Vv), 
the  result  is  the  thermal  energy  equation: 

p-jgg  =  a: Vv  -  V* q  (A.1 


Substitution  of  the  expression  for  a  given  earlier 


results  in  the  following  form  of  the  thermal  energy  equation: 

-pV*v-V*q  +  l>  (A.  19) 

A.  ( V  •  v )  2  +  2 y  ( e  n2  +  e^2  +  e  332  )  +  4u(ei22  +  E  232  +  e  3i2  ) 

where  the  e's  are  the  components  of  the  rate  of  strain  tensor 
(see  ref  1  and  10).  4>  represents  the  dissipation  of  energy 

due  to  viscous  effects;  in  many  practical  applications  it 
is  assumed  to  be  negligible. 

In  terms  of  enthalpy  (h  =  e  +  p/p),  and  when  viscous 
dissipation  is  negligible,  the  thermal  energy  equation 
becomes 


De 

PDt 

4>  = 


Dh  _  Dp 
pDt  Dt 


V  •  q 


(A. 20) 


For  a  perfect  gas,  Dh  =  c  DT.  If  q  is  linearly  depen¬ 
dent  on  the  gradient  of  temperature,  the  thermal  energy 
equation  takes  the  form 


DT  _  Dp  . 
pCp  Dt  Dt  + 


V* (kVT)  (q 


■  k  V  T ) 


(A.  21  ) 


For  the  case  of  constant  properties  and  pressure,  this 


reduces  to  the  familiar  form 


Appendix  B.  The  Thermal  Mixing  Region 

Consider  the  flow  situation  depicted  in  figure  B.l. 

A  fluid  is  moving  in  the  x-ccordinate  direction  with  velocity 
V  .  The  fluid  below  the  x-axis  is  initially  at  temper¬ 
ature  T  and  the  fluid  above  the  x-axis  is  held  at  T  . 

c  °° 

Diffusion  of  enthalpy  (or  temperature,  for  the  constant 
property  case)  is  prohibited  until  the  fluid  reaches  the 
position  x  =  0. 

At  this  point,  a  temperature  mixing  region  will  begin 
to  develop  and  grow  until  the  diffusion  (conduction)  of 
energy  is  balanced  by  convection.  This  balance  is  the  con¬ 
dition  for  steady-state. 

The  extent  of  the  temperature  mixing  region  must  be 
determined.  Clearly,  the  width  of  the  mixing  region,  2b, 
can  be  expected  to  grow  with  increased  distance  from  the 
point  where  mixing  begins.  As  a  first  cut  approximation, 
one  might  expect  isotherms  as  plotted  in  the  x-y  plane  to 
have  the  form  y  =  kx  (straight  lines).  As  the  mixing 
region  grows  in  the  downstream  direction,  however,  the 
heat  flux  at  an  isotherm  will  decrease.  This  suggests 
that  the  rate  of  growth  of  the  temperature  mixing  region 
should  decrease  with  increased  distance  from  the  initial 
point  of  mixing.  A  more  appropriate  relationship  for 

/  -i  /  "I 

isotherms  might  therefore  be  assumed  as  y  =  kx^ 

For  the  constant  property  case,  the  energy  equation 


is  equation  A. 22: 


which  is  expanded  in  terms  of  two-dimensional  Cartesiar 
coordinates  as 


3T  ,  3T  ,  3T  / 3 2  T  ,  32T- 

n  +  +  VW  =  (IP-  +  W 


The  following  assumptions  will  be  used  in  reducing 


equation  B.2: 


a)  steady-state  conditions  exist,  so  that  =  0. 

b)  The  momentum  (velocity)  field  is  unaffected  by 
the  temperature  mixing  region,  so  that  u  - 
and  v  =  0. 

c)  Gradients  in  the  y-direction  are  much  larger  than 

3  2  T 

gradients  in  the  x-direction  so  that  -g^r  is  neg- 

3  2  T 

ligible  compared  to  (the  "boundary  layer 

assumption" ) . 


Under  these  conditions,  equation  B.2  becomes 


3T  =  _a , 3 2  T  v 
3  x  V  [TyTl 

CO  v 


It  is  generally  convenient  to  non-dimensionali ze  the 
variable  of  interest  in  equations  such  as  B.3,  and  there¬ 
fore  a  non-dimensional  temperature  difference,  6,  is 


defined  as 


and  equation  B.3  becomes 

30.  _a_/  326  \ 

3x  V 

CO  V 

The  boundary  conditions  which  apply  to  the  steady- 
state  temperature  field  are: 

for  y  >  0,  T ( 0 , y )  = 
for  y  <  0,  T( 0  ,y)  =  Tc 
as  y  T(x,y)  -►  T^ 

as  y  +  -00,  T(  x , y )  •>  Tc 
T(x,0+)  =  T( x , 0" )  =  i(Tc  +  Tj 

Under  the  transformation  to  the  non-dimensional  var 
able  0  these  conditions  become: 

for  y  >  0,  0(0, y)  =  0 
for  y  <  0,  6(0, y)  =  1 
as  y  °°>  0  (  x , y )  0 

as  y  *  -oo,  9  ( x,y)  1 
0(x,O+)  =  0 ( x , 0” )  =  5 


The  problem  described  by  figure  B.1  is  therefore 
defined  by  equation  B.5  and  boundary  conditions  B.7  Two 


possible  methods  of  solution  are  convenient  for  this  prob¬ 
lem:  the  use  of  Laplace  transforms  on  the  x-variable  to 

obtain  an  ordinary  differential  equation  in  the  Laplace 
domain,  or  the  determination  of  a  similarity  transformation 
(x,y)  -*■  (q)  to  convert  equation  B.5  to  an  ordinary  differen 
tial  equation  in  q •  Both  of  these  methods  are  examined 
below . 

The  Laplace  Transform  Method 

Use  of  the  Laplace  transform  in  this  case  is  motiv¬ 
ated  by  the  fact  that  the  energy  equation  ■  1  -  La;  ! 

ian  ,  V2T.  The  existence  of  exhaust  i  >;•  •  .  -  :  '  !  •,  1  '  • 

tables  of  Laplace  transforms  and  inverse  Laplace  tran  s  s 

suggests  that  a  solution  of  the  transformed,  ordinary 
differential  equation  will  be  readily  invertible,  resul¬ 
ting  in  a  closed-form  solution  in  the  original  (x,y) 
space. 

Denoting  by  0  the  Laplace  transform  of  6  (and  trans¬ 
forming  on  the  variable  x) ,  then 

0(s,y)  =  j°°  e_sx9(x,y)  dx  (R. 

o 

Applying  the  transformation  to  equation  P.R  and 
boundary  conditions  B.7  yields  the  following  ordinary 
differential  equation  and  associated  boundary  conditions: 


for  y  >  0:  s0  =  — (jrr) 

CO  y 

as  y  0  ->■  0 

0(s,O+)  =  1/ ^  2s) 

for  y  <  0:  s0  =  +  1 

CO  V 

as  y  ,  0  ■>  1/s 

0 ( s , 0” )  =  1 /(2s) 


The  solution  for  the  case  y  >  0  must  be  established  sep¬ 
arately  from  the  solution  for  the  case  y  <  0  due  to  the 
fact  that  0(x,y)  has  different  initial  conditions  on  the 
variable  x  for  these  two  cases.  However,  the  condition 
0(s,O+)  =  0(s,O  )  assures  continuity  of  the  solution  for 
all  values  of  y. 

For  y  >  0,  the  solution  to  the  second  order,  linear 
differential  equation  has  the  general  form 


0  (  s  , y )  =  A  i  (  s)  exj 


Bi  ( s) exp( 


The  condition  for  convergence  as  y  +  »>  requires  that 
Ai(s)  =  0,  and  the  condition  on  0(s,O+)  requires  Bi(s)  = 
1/(2s).  Therefore, 

0  (  s  ,  y )  =  { 1  /  (  2s) }  exp( -y//s\F7a)  (y  >  0)  (B. 

For  the  case  y  <  0,  the  homegeneous  solution  is 


as  given  in  (B.10),  above,  and  a  particular  solution  is 
0  (s,y)  =  1/s.  The  general  solution  for  this  case  is  then 


S(s,y)  =  1/s  +  A 2 ( s) exp( y/ / sV^/ a )  +  B2 ( s) exp( -y/ZsV^/a) ( B. 1 2' 

Here,  the  condition  for  convergence  as  y  +  -»  requires 
B2(s)  =  0.  The  condition  on  8(s,0~)  requires  A2(s)  = 

-l/(2s).  Therefore, 

9 ( s ,y)  =  1/s  -  { 1 / ( 2s) } exp( - | y | / /sV^/a)  ( y  <  0)  ( B  .  1 3 ) 

The  inverse  Laplace  transforms  of  these  functions  are 
( ref  13): 

9(x,y)  =  ^erfc(y//4ax7T^)  (y  >  0)  (B.14) 

0 ( x , y )  =  1  -  2 erf c ( | y | //4ax/ V^)  (y  <  0)  ( B .  1 5 ) 

where  erfc(x)  =  1  -  erf(x)  and  erf(x)  is  given  by 

erf(x)  =  -j-j  Xexp(  -u2  )  dw  (B.16) 

v  IT  0 


Determination  of  a  Similarity  Transformation 


As  noted  earlier,  it  seems  reasonable  to  expect 
isotherms  to  be  plotted  in  the  (x,y)  plane  as  y  =  kx 


(1/n) 


In  such  a  case  it  is  possible  to  reduce  the  partial  differ¬ 
ential  equation  in  x  and  y  to  an  ordinary  differential 
equation  in  terms  of  a  similarity  transformation  (x,v)  ■+ 


p .  To  determine  the  appropriate  functional  relationship 
between  p  and  (x,y),  the  solution  form  is  assumed  to  be 


a  b/  a  s  c 
H  =  x  y  (  — ) 


(B 


Using  primes  (')  to  denote  differentiation  with  respect 
to  n>  we  have 


86  _  jj_0.  _9p.  _  ,  ar]_ 

9x  9p  3x  x 


90  _  90  9n  _  a  -  bp 

37  =  ~  ~  =  6  " 


9p  9y 


y 


(B 


9  2  0 


_  _90_  /  p,  |  bp  v  _9p 
9y2  9n  1  y  ;  3y 


-  0 


m  bln 


Substituting  these  expressions  into  equation  B.5  yield: 


nn9.  =  Wv^)b2n29„  +  n9 , (l,^./v?.)bi  _  „(a,/,v¥)b) 

x  y 


y 


y 


(B. 


Multiplying  this  expression  by  (— )  and  simplifying  yield: 

a. 


n0»  =  (^i)(xW/V<x.))n26„  +  (b)(bx(a/voo)  _  x(a/VcJ)ne, 

a  y  a  y  y 

( B , 


Inspection  of  equat ion  B.20  suggests  a  similarity 
transformation  of  the  form  px  =  y2/(ax/Voo),  in  which  case 
equation  B.5  reduces  to  the  ordinary  differential  equation 


4q  1 9 "  +  (2  +  pi)0'  =  0  (Pi 


y 2 / (ax/V  ) )  (B. 


There  is,  however,  a  somewhat  more  convenient  (and 
more  conventional)  choice  for  n •  If  b  is  chosen  equal  to 
unity,  then  the  last  term  on  the  right  hand  side  of  equation 
B.19  is  zero.  The  remaining  term  on  the  right  hand  side  of 
the  equation  must  have  (otx/V  )/y2  equal  to  for  self¬ 
similarity.  Since  t>  (the  power  of  y)  has  already  been 
chosen  to  be  one,  p  must  equal  -2. 

An  additional  simplification  occurs  if  n 2  is  chosen 
to  be  (ay)2/(ax/Vm)  (i.e.  n=  ± ( ay) / /ax/ Vot ,  with  the  factor 
(a)  absorbed  into  the  non-dimensional  variable).  By  the 
construction  of  q ,  a  and  c  must  be  set  equal  to  -h  so  that 
q  =  ±y//4otx/ Vot .  In  this  case,  equation  B.20  reduces  to 

0"  +  2r)0  '  =0  (n  =  ±y/ /4oix/ Vro)  (B.22) 

The  solution  to  equation  B.22  has  the  form  0  =  A  +  Berf(n). 
Application  of  the  boundary  conditions  (B.7)  again  yields 
equations  B . 1 4  and  B.15  as  the  required  solution  to  the 
differential  equation  B.5. 


Appendix  C.  Condensation  and  !_ 


The  equation  of  motion  for  a  particle  moving  with 
respect  to  a  fluid  may  be  obtained  by  determining  the- 
forces  which  act  on  the  particle  and  applying  Newton's 
Second  Law. 

In  general,  three  forces  will  act  on  such  a  parti  cl 
a  buoyant  force,  equal  in  magnitude  to  the  weight  of  trio 
displaced  fluid;  a  gravitational  force,  equal  ir.  magnitu 
to  the  weight  of  the  particle;  and  a  viscous  drag  force, 
which  depends  on  the  geometry  of  the  particle  under  con¬ 
sideration  and  is  directed  to  oppose  the  velocity  of  the 
particle  relative  to  the  fluid. 

If  the  vector  representing  the  acceleration  due  to 
gravity  is  assumed  to  be  parallel  to  the  y-axis,  these 
forces  can  be  modelled  by: 


Buoyant  force: 
Gravitational 
Viscous  Drag: 


force: 


x  Vol  x  g) j 

-(p  x  Vol  x  g)j  (i 


-(CD(ipfvr 


vr)A> 


V 

r 


where  the  subscripts  f  and  p  refer  to  i roper  lies  of  the 

fluid  and  the  particle ,  respectively;  Vol  is  the  volume 

of  the  particle  and  A  is  its  frontal  area;  v  _  is  the 

velocity  of  the  particle  relative  to  the  fluid  and  v  i 
1  r 

a  unit  vector  aliened  with  v  ;  and  1.,  is  the  nor:  -  i  :  <  e. - 

r 

sional  coefficient  of  drag  based  on  frontal  nr- -a  . 


The  value  of  is  dependent  on  the  geometry  of  the 
particle  and  on  the  Reynolds  number  of  the  particle  with 
respect  to  the  fluid. 

If  the  particle  can  be  approximated  as  a  sphere  then 
0^  as  a  function  of  Re  can  be  determined  from  the  plot  foun 
on  page  17  of  ref  1.  For  the  purpose  of  integrating  the 
equation  of  motion  for  the  particle  it  is  most  convenient 
to  establish  an  empirical  functional  relationship  between 
and  Re. 

Stoke' s  approximation  gives  =  24/Re  for  very  low 
Reynolds  numbers,  and  is  approximately  constant  for  the 
range  2000  <  Re  <  200,000  ( =  0.405).  In  the  interme¬ 
diate  range,  0.25  <  Re  <  2000,  the  data  suggests  a  poly¬ 
nomial  relationship  between  log^ q  and  log^Q  Re. 

Therefore,  relationships  of  the  form 


g  --  A54  +  B5 3  +  04 2  +  DC  +  E 
g  =  F52  +  G5  +  H 


5  =  log10  CD 


(0, 


1°*1 0  Re 


are  assumed,  and  a  least-squares  curve  can  be  fitted  t o 
the  data  presented  in  ref  1.  The  results  of  these  curvi 
f  i  t  s  a  r  e : 


n.0011 If4  -  0.1 405 3  +  0.1 150 2 
0.1 02g 2  -  0.8945  +  1.42 


0.88 8  *  4  1.42 


Table  C.1  compares  the  values  of  fro-:,  ro:'  1  with 
those  obtained  using  the  quadratic  and  quartic  polynomial 
approximations  given  above.  In  both  cases,  agreement  bet we 
the  accepted  value  and  the  approximation  is  very  rood  (with 
in  three  per  cent  in  most  cases)  . 

For  the  situation  under  consideration  here,  F.e  is  to 
be  evaluated  using  the  magnitude  of  v  .  With  this  relation 
ship  between  v  and  ,  the  equation  of  motion  for  a  sphere 
can  be  numerically  integrated  to  obtain  v  (t)  and  again  to 
determine  the  position  vector,  r(t). 

Let  the  velocity  of  the  particle  be  denoted  by  v(t) 

=  u(t)i  +  v(t)j  (note  that  v(t)  is  not  the  relative  velocit 
between  the  particle  and  the  fluid) .  If  the  spherical  par¬ 
ticle  is  assumed  uo  have  been  formed  instantaneously  at  tin 
t  =  0,  then  v(0)  =  V  ,  where  is  the  vector  representing 
the  fluid  velocity.  For  the  case  where  the  free  stream 
velocity  is  perpendicular  to  the  gravitational  acceleration 
vector,  V  =  V  i . 

CO  00 

A  qualitative  argument  establishes  the  fact  that 
u(t)  =  V  under  these  conditions.  Since  the  drag  force 

'X' 

depends  on  the  relative  velocity  between  the  particle  and 
the  fluid,  there  can  be  no  viscous  drag  on  the  particle  at 
‘  t  =  0.  The  gravitational  and  buoyant  forces  do  rot 

:'f  oei  the  x-oonoonent  of  velocity,  so  u(  At)  must  equal 
At  time  t  =  At,  therefore,  v  must  be  nor non i i c - 
*  ■  the  free  stream  velocity  with  no  component  '  the 
This  argumen  t  can  be  extended  to  any  t  i 


71 


t  to  show  that,  for  V  =  V  1,  there  can  be  no  x-eo.\  r  on-nt 

CO  oo 

to  the  vector  v  . 

r 

Furthermore,  this  yields  the  result  that  v  =  ± j ,  and 
since  the  spherical  particle  is  expected  to  fall  under  the 
influence  of  gravity,  v^  =  -j. 

Therefore,  for  the  case  under  consideration  the 
equations  of  motion  are  given  by 

u  ( t )  =  V 

p  x  Vol  x  v  =  (p  -  p  )  x  Vol  x  g 

p  f  p  u 

+  CD(sPfv2A)  (C 

v( 0)  =  0 

The  volume  of  a  sphere  is  given  by  Varrr3  ,  and  its 
frontal  area  is  vr2 .  Using  these  relationships  in  the 
equation  above  and  dividing  both  sides  by  p^  x  Vol  yields 
the  following  result: 

u  ( t )  =  V 

v  =  (p-l)g  +  GD(3/8r)pv2  (C 

v  (  0 )  =  0 

where  p  is  defined  by  p  =  p^/p^.  It  is  also  convenient 
to  define  the  parameter  f"  to  be  g*  =  (p-l)g.  A  good 
approximation  to  this  value  is  g"  -  32.1  ft/sec2. 

When  the  acceleration  dun  to  effective  gravitational 
force  ( g ” )  is  balanced  by  the  viscous  drag  force’,  the 


spherical  particle  will  fall  at  a  constant  velocity  for 
all  time  from  that  point  on.  This  possibility  must  be 
accounted  for  before  attempting  a  numerical  solution  to 
equation  C.6.  For  this  purpose,  the  condition  v  -  0  if 
G^(3/8r)pv2  >  g*  is  added  to  equation  0.6. 

Several  types  of  numerical  schemes  were  used  to 
integrate  equation  C.6.  /.'her.  these  schemes  wer c  compared 

it  was  determined  that  a  simple  forward  differencing  tech¬ 
nique  gave  sufficient  accuracy  in  its  results.  The  for¬ 
ward  difference  equations  employed  in  this  algorithm  are: 

v ( t+At)  -  v ( t )  +  At(CD(3/8r)pv2  -  32.1),  if 

CD(3/8r)pv2<  32.1 

v(t+At)  =  0,  otherwise  (C.7) 

v(0)  =  0 


For  the  purposes  of  this  study,  the  fluid  properties 
were  taken  to  be  those  of  nitrogen  at  -280  F  and  1  atm 


pressure  and  the  value  of  p  was  chosen  to  be  the  density 
of  water  at  32  F. 


The  numerical  integration  was  performed  on  a  Hewlett- 
Packard  model  9845  desktop  computer  using  a  program  written 
in  BAfIC.  The  results  of  tho  numerical  integration  for 


some  values  of  r  and  At  are  presented  in  tho  tables  which 


follow. 


Tab 
cal 
t  = 


| 


S3 


t  ( sec ) 

Distance  Travelled  (ft) 

2r=0 . 05  ft 

2r=0 . 01  ft 

2r=0.005  ft 

0.1 

0.16 

0.16 

0.16 

0.2 

0.64 

0.64 

0.64 

0.3 

1.4 

1.4 

1.4 

0.4 

2.5 

2.4 

2.3 

0.5 

3.9 

3.7 

3.4 

0.6 

5.6 

5.1 

4.5 

0.7 

7.6 

6.6 

5.7 

0.8 

9.7 

8.2 

6.9 

0.9 

12 

9.8 

8.1 

1  .0 

15 

1 1 

9.4 

Table  C .  3 . 

Results  of  the  Numerical  Integration  of  Equation  C.7 

At  =  0.01  sec 


t  (  sec ) 

Distance  Travelled  (ft) 

2r =0 . 05  ft 

2r=0.01  ft 

0.01 

0.0016 

0.0016 

0.02 

0.0064 

0.0064 

0.03 

0.014 

0.014 

0.04 

0.026 

0.026 

0.05 

0.040 

0.040 

0.06 

0.058 

0.058 

0.07 

0.079 

0.078 

0.08 

0.10 

0.10 

0.09 

0.1  3 

0.13 

0.10 

0.16 

0.16 

Appendix  D.  Property  Values  and  Sample  Calculations 

In  order  to  determine  the  values  for  the  various  par¬ 
ameters  which  have  been  shown  to  be  important  in  this  Rey¬ 
nolds  number  simulation  technique  (e.g..  Re,  L) ,  it  is 
necessary  to  obtain  various  property  values  as  functions 
of  pressure  and  temperature.  In  many  cases,  these  values 
can  be  read  from  tables  (see,  for  example  refs  14  and  15). 
For  convenience,  expressions  which  can  be  used  to  evaluate 
many  of  these  properties  are  summarized  below. 


Density 


From  the  ideal  gas  law,  the  value  for  density  can 
be  determined  from  p=  -j^,  where  R  is  the  gas  constant  for 


the  particular  gas. 

The  value  of  R  is  determined  from  R  =  R/M,  where  R 
is  the  universal  gas  constant  (1545.43  ( f t-lb^/R-lb^  mole)) 
and  M  is  the  molecular  weight  of  the  gas.  For  air,  R  - 
1716  ft2/R-sec2  and  for  nitrogen,  R  -  1775  ft2/R-sec2. 


Viscosity 


The  coefficient  of  viscosity,  u»  is  not  strongly 
dependent  on  pressure.  An  empirical  relationship  between 


U  and  T  is  given  by  Sutherland's  law  to  within  ±2!  when 
T  is  within  prescribed  limits.  The  relationship  is 


U.  ( X  \/2  T  o  t  S 
Uo  T0  ^  T  +  S 


(  D .  1  ) 


For  nitrogen:  Uo  =  3.473x10-7  slug/ft-sec,  T0  =  491.6  R 
S  -  192  R,  and  180  R  <  T  <  2700  R. 

For  air:  Uo  =  3.584  x  10~7  slug/ft-sec,  T0  =  491.6  R, 

S  =  199  R,  and  300  R  <  T  <  3420  R. 

The  value  for  the  kinematic  viscosity,  v,  can  be 
determined  from  the  relationship  v  =  u/p. 


Thermal  Conductivity 


A  Sutherland's  law  relationship  also  exists  for  the 
coefficient  of  thermal  conductivity,  k. 


k  ,T  J/2  T0  +  S 
k  o  T0  T  +  S 


(D.2) 


For  nitrogen:  k0  =  0 . 01 40 B/h-ft-R,  T0  =  491.6  R, 

S  =  300  R,  and  260  R  <  T  <  2160  R. 

For  air:  k0  0 . 01  395  B/h-ft-R ,  T0  =  491.6  R, 

5  =  350  R,  and  300  R  <  T  <  1800  R. 

The  value  for  thermal  diffusivity,  a,  can  be  derived 
from  nt  =  k/p c  . 


Sample  Calculation 

In  order  to  illustrate  the  use  of  the  techniques  out 
lined  in  the  paper  to  simulate  high  Reynolds  number  flow, 
a  sample  calculation  is  carried  out  below: 

It  is  desired  to  simulate  Re  ~  IxlO7  on  a  flat  plate 

two  feet  long  and  one  foot  wide.  The  tunnel  total  temper 

ature  is  approximately  540  R,  and  total  pressure  (for  the 

tunnel  air)  is  about  1  atm  (2117  lb^/ft2). 

Liquid  nitrogen  is  used  as  the  cool  fluid,  so  its 

properties  are  evaluated  at  =  180  R.  From  ref  14, 

uc  =  4.6l1xl0-6  lb^/ sec -ft  =  1.432x10'7  slug/sec-ft 

k  =  5.460x10'8  B/h-ft-R 
c 

At  a  pressure  of  1  atm,  the  density  of  nitrogen  is 

0.2173  lb  /ft3;  this  value  can  be  used  to  determine  the 
m 

required  V^. 

Re  =  p  V  £/u  =  IxlO7,  so 

Mc  <»  c  ^ 

Vro  =  (lx107)/(pc£/uc) 

The  quantity  pc-£/uc  is  approximately  9.42x104,  so 
the  required  velocity  is  V  -  1 06  ft/sec. 

Chapter  V  outlines  the  procedure  for  determining  p 
from  ,  T^,  and  p^ . 

The  density  of  the  tunnel  air,  p  ,  is  assumed  equal 


2 


q  =  ip  V 

n  s  00 


-  q 


H 2 . 3x1 0~ 3  ) ( 1 06 )  2  =  12.92  lbf/ft 
2117  -  12.92  *  2104  lbf/ft 


ps/pt  =  0*994 


(D.4 


For  air  at  low  temperatures,  y  =  1.4*  and  the  ratio 
( Y  —  1 ) / y  -  0.286.  The  isentropic  relations  give 


Ts  =  0.998  Tt  =  T  =  540  R 


(0.5 


and  y'(Ts  =  540)  =  y  -  1.4,  so  Tg  -  540  R. 

The  value  of  p  '  is  found  from  the  ideal  gas  law. 
s 

Pc'  =  p  /K  T  =  2104/(1716  X  540)  =  2 . 27x 1 0_  3  slugs/ft 

s  s  s  s 

-  ps  (0.6 

Since  p  =  p  ' ,  p  can  be  taken  as  2104  lb~/ft2,  and 

S  S  S  X 

therefore  p  =  p  =  2104  lb„/ft2. 

C  S  1 

It  is  now  possible  to  determine  the  values  vc  and  a  . 

Pc  =  pc/RcTc  =  2104/(1775  x  180)  =  6.58x10“3  slugs/ft3 
v  =  y  /p  =  1 .432x10“7/6.58x10'3=  2 . 1 76x 1 0~ 5  ft2/sec 

C  C  c 

a  =  k  /p  c  =  5*46x1 0‘3/(6.58xi0"3x0. 255) 
c  cep 

F  c 

-  3*254  lb  -ft2 / slug-hr  -  0.101  ft2/hr 
m 

=  2 . 807x 1 0- 5  f t2 / sec 


where  the  value  of  c  has  been  taken  as  0.255  B/lb  -R. 


With  these  values  for  a  and  v  ,  it  is  possible  to 
determine  the  value  of  the  parameter  L.  Since  the  boundar 
layer  is  expected  to  be  turbulent  for  Re  =  1><107,  the 

_____  _____  3/ 

value  to  be  calculated  is  L.  =  3-3/a  / V  +  0.37/v  /V  l'10. 

£  c  00  c  00 

/2. 8  07x1  0  5  /I  06  =  5- 146x1  0"  4  ft^2 

?2. 176x10  5  / 1  06  =  4- 597x10" 2  ft^5  (D 

3/  3/ 

2/l0=  1.23  ft/l0 

1/ 

This  gives  L ^  -  2.262x10  2  ft/2.  The  constraint  on 
h  requires  h  >  and  -  0.032  ft,  so 

h  >  0.032  ft  (D 


/a  /V 
c  ° 


/V 

C  0 


•7l0  _ 


To  determine  the  required  mass  rate  of  flow,  the 
applicable  constraint  is  m/w  >  Pc7ooi^.i/Z,  so 

m/w  >  6.58x10“3  x  106  x  0.032  =  2.23><10"2slugs/ft-sec 

(D 

For  the  plate  under  consideration,  w  =  1  ft,  so 


Ik  >  2.23x1  0  2  slurs/sec  -  0.719  lb  /sec 

m 


If  condensation  or  icing  is  expected  to  be 
the  value  t/V  should  be  calculated  and  then  us 

OO 

if  condensation  will  effect  the  boundary  layer. 
For  this  exercise,  t  /  V  -j  0.0  2  s  o  c  . 

'  CO 


a  problem, 
e  d  t  o  s  e  e 


oa  r  - 


A  s  pii  e  r  i  c  a  1 


tide  with  a  diameter  of  0.01  ft  (  about  %  in)  will  fall 

approximately  0.006  ft  (0.07  in)  in  this  time.  For  the  ease- 
under  consideration,  k  >>  0.006  ft,  so  condensation  and  Icin 
are  not  expected  to  effect  the  Reynolds  simulation. 

The  value  for  the  Reynolds  number  can  be  determined  fro 
the  flow  properties  which  have  been  established  above. 

Rec  =  Voo£/vc  = ( 1 06) ( 2) /2 . 1  76x1  0” 5  -  9.74*106  (D.11 

This  is  very  close  to  the  desired  value  of  fie  =  1*107. 

It  is  instructive  to  calculate  the  value  of  Re  ,  the 

s 

Reynolds  number  based  on  properties  of  the  tunnel  air.  It 

has  already  been  established  thatpr,  =  2.27x10  3slugs/ft3, 

and  u  is  found  to  be  3.834x10-7  slugs/sec-ft .  Therefore, 
s 

v  -  1 .698x10" 4  ft2 /sec,  and 

Re„  =  Va£/v^  =  (106) (2)/1 .698xl0'4  =  1.25*106  (D.12 

From  this,  the  ratio  P.e  /Reo  can  be  found  to  be  7.8. 

The  cooling  technique  described  in  this  paper  has  therefore 
roulted  in  an  increase  in  effective  Reynolds  number  of  nearl 
an  order  of  magnitude. 


A  p  p  e  n  d  i  x  E  . 


The  Parameter  L 


The  parameter  L  which  appears  in  the  statement  of 
constraints  developed  in  Chapter  V  contains  all  of  the  inforn 
ation  regarding  the  temperature  dependency  of  those  con¬ 
straints  . 

This  temperature  dependency  can  be  made  explicit 
by  applying  the  Sutherland's  law  equations  presented  in 
Appendix  D  along  with  the  ideal  gas  law . 

From  the  ideal  gas  law, 


£ 

0  o 


(*)(£*> 

l  p  o 


<¥><£> 


(E.1  ) 


where  T0  is  taken  as  491-6  R,  p0  =  1  atm,  and  p0  =  p0/RT0. 
If  v o  is  defined  to  be  uo/po>  then  Sutherland's  law 

g  i  v  e  s 


v 

Vo 


{  (y~)5/2 
1  o 


T 


±JL 


■f  <f> 


(E.2) 


A  similar  expression  can  be  developed  for  the  thermal 
diffusivity,  a.  If  c  is  assumed  constant,  and  for  o0  - 
ko/PoC  ,  >3  result  is 


a 

a0 


T 


La  +  £.1)  (iia.) 

T  +  s*'  k  p  ' 


3) 


The  values  for  S,  f>“ ,  and  T0  are  given  in  Appendix  D 


In  the  case  where  the  tunnel  pressure,  t 
the  same  as  atmospheric  pressure,  the  quantit 


nearly  unity.  With  this  assumption,  and  fro::, 
of  L,  the  temperature  dependence  of  L  is  expr 


(i)s/'(/S7/r;( 
±  0 


To  +  S--sV; 
°°XT  +S 


)/2  } 


1/ 


(^)'2{^77v^  +; 


+  f 


)'Sli 


From  equations  E.4  and  E.5  it  is  apparer. 
an  increasing  function  of  temperature,  as  wel 
ing  function  of  free  stream  tunnel  velocity, 
teristics  are  illustrated  in  the  figures  that 
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